set.seed(1234)
library(splines)
library(ggplot2)
library(viridis)
library(patchwork)
library(hexbin)
library(reshape2)
library(ResistPhy)
library(ape)
library(posterior)
library(cmdstanr)
library(bayesplot)
library(latex2exp)

run_mcmc <- T

Ingest synthetic usage

We Restrict ourselves to a 20 year span Load usage data and simulation script

usage_df <- simulated_usage()
source("simu_phylo.R")

Set simulation parameters and simulate phylogenies

q_u <- 1.2
q_t <- -2.7

n_tip <- 200

gamma_sus <- 1/60.0*365.0

gamma_res_u <- gamma_sus + q_u
gamma_res_t <- gamma_sus + q_t

t0 <- min(usage_df$time)
tmax <- max(usage_df$time)
sim <- sample_phylo(gamma_sus, gamma_res_u, gamma_res_t, n_tip, 123456, 1234567) 

df <- sim$traj
phys <- sim$phys
mr <- sim$most_recent_samp

Visualise (t)

N_0 <- 3e6
incidence <- 1.5e5 / N_0
mu <- incidence+gamma_sus

time_grid <- seq(from=t0, to=tmax, by=1e-4)
t_change <- (t0 + floor(0.33*(tmax-t0)))

beta_func  <- function (t) if (t < t_change) 1.02*mu else mu*(0.98 + 0.06*(t-t_change)/(tmax-t_change))
N_func <- function(t) N_0*exp(log(2)*(t-t0)/(tmax-t0))
usg_func <- yearly_usg_stepfunc(usage_df$usage, usage_df$time)

beta_vals <- sapply(time_grid, beta_func)
beta_df <- data.frame(time=time_grid, beta=beta_vals)

usg_vals <- sapply(time_grid, usg_func)
usg_df <- data.frame(time=time_grid, usg=usg_vals)

N_vals <- sapply(time_grid, N_func)
N_df <- data.frame(time=time_grid, N=N_vals)

N_plt <- ggplot(N_df, aes(x=time, y=N))+
    geom_line() +
    xlim(t0,tmax) +
    labs(y=TeX(paste0("$","N(t)","$")))+
    theme_minimal()+
    theme(axis.title.x=element_blank(),
        axis.text.x=element_blank(),
        axis.ticks.x=element_blank())
beta_plt <- ggplot(beta_df, aes(x=time, y=beta))+
    geom_line()+
    xlim(t0,tmax) +
    labs(y=TeX(paste0("$","\\beta(t)","$")))+
    theme_minimal()+
    theme(axis.title.x=element_blank(),
        axis.text.x=element_blank(),
        axis.ticks.x=element_blank())
usg_plt <- ggplot(usg_df, aes(x=time, y=usg))+
    geom_line()+
    xlim(t0,tmax) +
    labs(x="Time",y=TeX(paste0("$","u(t)","$")))+
    theme_minimal()

f_plt <- N_plt/beta_plt/usg_plt
pdf("Figures/functions.pdf",4,6)
f_plt
dev.off()
#> png 
#>   2
f_plt

Visualise usage

usg_fun <- yearly_usg_stepfunc(usage_df$usage, usage_df$time)
ggplot(data.frame(time=time_grid, usage=sapply(time_grid, usg_fun)), aes(x=time, y=usage))+geom_line()+theme_minimal()

plot(ggplot(df, aes(x=t))+geom_line(aes(y=I_2),color="blue") +geom_line(aes(y=I_3),color="red"))

Preview the phylogenies run BNPR to see how that looks

plot(phys[[1]])
axisPhylo()

plot(phys[[2]])
axisPhylo()

Run MCMC

out <- infer_costs2(phys, 
                    mr,
                    usage_df$usage,
                    usage_df$time,
                    t0,
                    tmax,
                    1/60.0,
                    365,
                    n_iter=2000, 
                    n_warmup=2000,
                    model="nodecay",
                    K=60,
                    L=6.5,
                    gamma_log_sd=0.3,
                    stan_control=list(adapt_delta=.99,
                        max_treedepth=13,
                        parallel_chains=4,
                        chains=4,
                        refresh=10
                    ))                    
#> Compiling Stan program...
#> 
-

\

|

/

-

\

|

/

-

\

|

/

-

\

|

/

-

\

|

/

-

\

|

/

-

\

|

/

-

\

|

/

-

\

|

/

-

\

|

/

-

\

|

/

-

\

|

/

-

\

|

/

-

\

|

/

-

\

|

/

-

\

|

/

-

\

|

/

-

\

|

/

-

\

|

/

-

\

|

/

-

\

|

/

-

\

|

/

-

\

|

/

-

\

|

/

-

\

|

/

-

\

|

/

-

\

|

/

-

\

|

/

-

\

|

/

-

\

|

/

-

\

|

/

-

\

|

/

-

\

|

/

-

\

|

/

-

\

|

/

-

\

|

/

-

\

|

/

-

\

|

/

-

\

|

/

-

\

|

/

-

\

|

/

-

\

|

/

-

\

|

/

-

\

|

/

-

\

|

/

-

\

|

/

-

\

|

/

-

\

|

/

-

 
Running MCMC with 4 parallel chains...
#> 
#> Chain 1 Iteration:    1 / 4000 [  0%]  (Warmup)
#> Chain 1 Informational Message: The current Metropolis proposal is about to be rejected because of the following issue:
#> Chain 1 Exception: log1p: x is -2.39349e+19, but must be greater than or equal to -1.000000 (in '/tmp/Rtmpb97qYm/model-94ed8f713.stan', line 197, column 8 to column 82)
#> Chain 1 If this warning occurs sporadically, such as for highly constrained variable types like covariance matrices, then the sampler is fine,
#> Chain 1 but if this warning occurs often then your model may be either severely ill-conditioned or misspecified.
#> Chain 1
#> Chain 1 Informational Message: The current Metropolis proposal is about to be rejected because of the following issue:
#> Chain 1 Exception: log1p: x is -2.48809e+19, but must be greater than or equal to -1.000000 (in '/tmp/Rtmpb97qYm/model-94ed8f713.stan', line 197, column 8 to column 82)
#> Chain 1 If this warning occurs sporadically, such as for highly constrained variable types like covariance matrices, then the sampler is fine,
#> Chain 1 but if this warning occurs often then your model may be either severely ill-conditioned or misspecified.
#> Chain 1
#> Chain 1 Informational Message: The current Metropolis proposal is about to be rejected because of the following issue:
#> Chain 1 Exception: log1p: x is -775.108, but must be greater than or equal to -1.000000 (in '/tmp/Rtmpb97qYm/model-94ed8f713.stan', line 197, column 8 to column 82)
#> Chain 1 If this warning occurs sporadically, such as for highly constrained variable types like covariance matrices, then the sampler is fine,
#> Chain 1 but if this warning occurs often then your model may be either severely ill-conditioned or misspecified.
#> Chain 1
#> Chain 1 Informational Message: The current Metropolis proposal is about to be rejected because of the following issue:
#> Chain 1 Exception: log1p: x is -1.01558, but must be greater than or equal to -1.000000 (in '/tmp/Rtmpb97qYm/model-94ed8f713.stan', line 197, column 8 to column 82)
#> Chain 1 If this warning occurs sporadically, such as for highly constrained variable types like covariance matrices, then the sampler is fine,
#> Chain 1 but if this warning occurs often then your model may be either severely ill-conditioned or misspecified.
#> Chain 1
#> Chain 2 Iteration:    1 / 4000 [  0%]  (Warmup)
#> Chain 2 Informational Message: The current Metropolis proposal is about to be rejected because of the following issue:
#> Chain 2 Exception: gamma_lpdf: Random variable is 0, but must be positive finite! (in '/tmp/Rtmpb97qYm/model-94ed8f713.stan', line 232, column 4 to column 24)
#> Chain 2 If this warning occurs sporadically, such as for highly constrained variable types like covariance matrices, then the sampler is fine,
#> Chain 2 but if this warning occurs often then your model may be either severely ill-conditioned or misspecified.
#> Chain 2
#> Chain 3 Iteration:    1 / 4000 [  0%]  (Warmup)
#> Chain 3 Informational Message: The current Metropolis proposal is about to be rejected because of the following issue:
#> Chain 3 Exception: gamma_lpdf: Random variable is 0, but must be positive finite! (in '/tmp/Rtmpb97qYm/model-94ed8f713.stan', line 232, column 4 to column 24)
#> Chain 3 If this warning occurs sporadically, such as for highly constrained variable types like covariance matrices, then the sampler is fine,
#> Chain 3 but if this warning occurs often then your model may be either severely ill-conditioned or misspecified.
#> Chain 3
#> Chain 3 Informational Message: The current Metropolis proposal is about to be rejected because of the following issue:
#> Chain 3 Exception: gamma_lpdf: Random variable is 0, but must be positive finite! (in '/tmp/Rtmpb97qYm/model-94ed8f713.stan', line 232, column 4 to column 24)
#> Chain 3 If this warning occurs sporadically, such as for highly constrained variable types like covariance matrices, then the sampler is fine,
#> Chain 3 but if this warning occurs often then your model may be either severely ill-conditioned or misspecified.
#> Chain 3
#> Chain 3 Informational Message: The current Metropolis proposal is about to be rejected because of the following issue:
#> Chain 3 Exception: gamma_lpdf: Random variable is 0, but must be positive finite! (in '/tmp/Rtmpb97qYm/model-94ed8f713.stan', line 232, column 4 to column 24)
#> Chain 3 If this warning occurs sporadically, such as for highly constrained variable types like covariance matrices, then the sampler is fine,
#> Chain 3 but if this warning occurs often then your model may be either severely ill-conditioned or misspecified.
#> Chain 3
#> Chain 3 Informational Message: The current Metropolis proposal is about to be rejected because of the following issue:
#> Chain 3 Exception: gamma_lpdf: Random variable is 0, but must be positive finite! (in '/tmp/Rtmpb97qYm/model-94ed8f713.stan', line 232, column 4 to column 24)
#> Chain 3 If this warning occurs sporadically, such as for highly constrained variable types like covariance matrices, then the sampler is fine,
#> Chain 3 but if this warning occurs often then your model may be either severely ill-conditioned or misspecified.
#> Chain 3
#> Chain 3 Informational Message: The current Metropolis proposal is about to be rejected because of the following issue:
#> Chain 3 Exception: gamma_lpdf: Random variable is 0, but must be positive finite! (in '/tmp/Rtmpb97qYm/model-94ed8f713.stan', line 232, column 4 to column 24)
#> Chain 3 If this warning occurs sporadically, such as for highly constrained variable types like covariance matrices, then the sampler is fine,
#> Chain 3 but if this warning occurs often then your model may be either severely ill-conditioned or misspecified.
#> Chain 3
#> Chain 4 Iteration:    1 / 4000 [  0%]  (Warmup)
#> Chain 4 Informational Message: The current Metropolis proposal is about to be rejected because of the following issue:
#> Chain 4 Exception: gamma_lpdf: Random variable is 0, but must be positive finite! (in '/tmp/Rtmpb97qYm/model-94ed8f713.stan', line 232, column 4 to column 24)
#> Chain 4 If this warning occurs sporadically, such as for highly constrained variable types like covariance matrices, then the sampler is fine,
#> Chain 4 but if this warning occurs often then your model may be either severely ill-conditioned or misspecified.
#> Chain 4
#> Chain 3 Iteration:   10 / 4000 [  0%]  (Warmup) 
#> Chain 2 Iteration:   10 / 4000 [  0%]  (Warmup) 
#> Chain 4 Iteration:   10 / 4000 [  0%]  (Warmup) 
#> Chain 2 Iteration:   20 / 4000 [  0%]  (Warmup) 
#> Chain 1 Iteration:   10 / 4000 [  0%]  (Warmup) 
#> Chain 4 Iteration:   20 / 4000 [  0%]  (Warmup) 
#> Chain 2 Iteration:   30 / 4000 [  0%]  (Warmup) 
#> Chain 4 Iteration:   30 / 4000 [  0%]  (Warmup) 
#> Chain 2 Iteration:   40 / 4000 [  1%]  (Warmup) 
#> Chain 4 Iteration:   40 / 4000 [  1%]  (Warmup) 
#> Chain 2 Iteration:   50 / 4000 [  1%]  (Warmup) 
#> Chain 1 Iteration:   20 / 4000 [  0%]  (Warmup) 
#> Chain 4 Iteration:   50 / 4000 [  1%]  (Warmup) 
#> Chain 2 Iteration:   60 / 4000 [  1%]  (Warmup) 
#> Chain 4 Iteration:   60 / 4000 [  1%]  (Warmup) 
#> Chain 2 Iteration:   70 / 4000 [  1%]  (Warmup) 
#> Chain 1 Iteration:   30 / 4000 [  0%]  (Warmup) 
#> Chain 4 Iteration:   70 / 4000 [  1%]  (Warmup) 
#> Chain 2 Iteration:   80 / 4000 [  2%]  (Warmup) 
#> Chain 1 Iteration:   40 / 4000 [  1%]  (Warmup) 
#> Chain 4 Iteration:   80 / 4000 [  2%]  (Warmup) 
#> Chain 2 Iteration:   90 / 4000 [  2%]  (Warmup) 
#> Chain 1 Iteration:   50 / 4000 [  1%]  (Warmup) 
#> Chain 4 Iteration:   90 / 4000 [  2%]  (Warmup) 
#> Chain 2 Iteration:  100 / 4000 [  2%]  (Warmup) 
#> Chain 1 Iteration:   60 / 4000 [  1%]  (Warmup) 
#> Chain 4 Iteration:  100 / 4000 [  2%]  (Warmup) 
#> Chain 3 Iteration:   20 / 4000 [  0%]  (Warmup) 
#> Chain 1 Iteration:   70 / 4000 [  1%]  (Warmup) 
#> Chain 1 Iteration:   80 / 4000 [  2%]  (Warmup) 
#> Chain 2 Iteration:  110 / 4000 [  2%]  (Warmup) 
#> Chain 1 Iteration:   90 / 4000 [  2%]  (Warmup) 
#> Chain 4 Iteration:  110 / 4000 [  2%]  (Warmup) 
#> Chain 1 Iteration:  100 / 4000 [  2%]  (Warmup) 
#> Chain 1 Iteration:  110 / 4000 [  2%]  (Warmup) 
#> Chain 3 Iteration:   30 / 4000 [  0%]  (Warmup) 
#> Chain 2 Iteration:  120 / 4000 [  3%]  (Warmup) 
#> Chain 1 Iteration:  120 / 4000 [  3%]  (Warmup) 
#> Chain 1 Iteration:  130 / 4000 [  3%]  (Warmup) 
#> Chain 2 Iteration:  130 / 4000 [  3%]  (Warmup) 
#> Chain 3 Iteration:   40 / 4000 [  1%]  (Warmup) 
#> Chain 4 Iteration:  120 / 4000 [  3%]  (Warmup) 
#> Chain 1 Iteration:  140 / 4000 [  3%]  (Warmup) 
#> Chain 3 Iteration:   50 / 4000 [  1%]  (Warmup) 
#> Chain 2 Iteration:  140 / 4000 [  3%]  (Warmup) 
#> Chain 1 Iteration:  150 / 4000 [  3%]  (Warmup) 
#> Chain 3 Iteration:   60 / 4000 [  1%]  (Warmup) 
#> Chain 2 Iteration:  150 / 4000 [  3%]  (Warmup) 
#> Chain 1 Iteration:  160 / 4000 [  4%]  (Warmup) 
#> Chain 3 Iteration:   70 / 4000 [  1%]  (Warmup) 
#> Chain 4 Iteration:  130 / 4000 [  3%]  (Warmup) 
#> Chain 2 Iteration:  160 / 4000 [  4%]  (Warmup) 
#> Chain 1 Iteration:  170 / 4000 [  4%]  (Warmup) 
#> Chain 3 Iteration:   80 / 4000 [  2%]  (Warmup) 
#> Chain 1 Iteration:  180 / 4000 [  4%]  (Warmup) 
#> Chain 3 Iteration:   90 / 4000 [  2%]  (Warmup) 
#> Chain 4 Iteration:  140 / 4000 [  3%]  (Warmup) 
#> Chain 2 Iteration:  170 / 4000 [  4%]  (Warmup) 
#> Chain 1 Iteration:  190 / 4000 [  4%]  (Warmup) 
#> Chain 3 Iteration:  100 / 4000 [  2%]  (Warmup) 
#> Chain 2 Iteration:  180 / 4000 [  4%]  (Warmup) 
#> Chain 1 Iteration:  200 / 4000 [  5%]  (Warmup) 
#> Chain 4 Iteration:  150 / 4000 [  3%]  (Warmup) 
#> Chain 2 Iteration:  190 / 4000 [  4%]  (Warmup) 
#> Chain 1 Iteration:  210 / 4000 [  5%]  (Warmup) 
#> Chain 2 Iteration:  200 / 4000 [  5%]  (Warmup) 
#> Chain 1 Iteration:  220 / 4000 [  5%]  (Warmup) 
#> Chain 3 Iteration:  110 / 4000 [  2%]  (Warmup) 
#> Chain 2 Iteration:  210 / 4000 [  5%]  (Warmup) 
#> Chain 1 Iteration:  230 / 4000 [  5%]  (Warmup) 
#> Chain 2 Iteration:  220 / 4000 [  5%]  (Warmup) 
#> Chain 1 Iteration:  240 / 4000 [  6%]  (Warmup) 
#> Chain 4 Iteration:  160 / 4000 [  4%]  (Warmup) 
#> Chain 2 Iteration:  230 / 4000 [  5%]  (Warmup) 
#> Chain 1 Iteration:  250 / 4000 [  6%]  (Warmup) 
#> Chain 2 Iteration:  240 / 4000 [  6%]  (Warmup) 
#> Chain 2 Iteration:  250 / 4000 [  6%]  (Warmup) 
#> Chain 3 Iteration:  120 / 4000 [  3%]  (Warmup) 
#> Chain 1 Iteration:  260 / 4000 [  6%]  (Warmup) 
#> Chain 3 Iteration:  130 / 4000 [  3%]  (Warmup) 
#> Chain 3 Iteration:  140 / 4000 [  3%]  (Warmup) 
#> Chain 4 Iteration:  170 / 4000 [  4%]  (Warmup) 
#> Chain 2 Iteration:  260 / 4000 [  6%]  (Warmup) 
#> Chain 1 Iteration:  270 / 4000 [  6%]  (Warmup) 
#> Chain 3 Iteration:  150 / 4000 [  3%]  (Warmup) 
#> Chain 1 Iteration:  280 / 4000 [  7%]  (Warmup) 
#> Chain 4 Iteration:  180 / 4000 [  4%]  (Warmup) 
#> Chain 1 Iteration:  290 / 4000 [  7%]  (Warmup) 
#> Chain 3 Iteration:  160 / 4000 [  4%]  (Warmup) 
#> Chain 4 Iteration:  190 / 4000 [  4%]  (Warmup) 
#> Chain 4 Iteration:  200 / 4000 [  5%]  (Warmup) 
#> Chain 1 Iteration:  300 / 4000 [  7%]  (Warmup) 
#> Chain 4 Iteration:  210 / 4000 [  5%]  (Warmup) 
#> Chain 1 Iteration:  310 / 4000 [  7%]  (Warmup) 
#> Chain 2 Iteration:  270 / 4000 [  6%]  (Warmup) 
#> Chain 4 Iteration:  220 / 4000 [  5%]  (Warmup) 
#> Chain 1 Iteration:  320 / 4000 [  8%]  (Warmup) 
#> Chain 4 Iteration:  230 / 4000 [  5%]  (Warmup) 
#> Chain 1 Iteration:  330 / 4000 [  8%]  (Warmup) 
#> Chain 4 Iteration:  240 / 4000 [  6%]  (Warmup) 
#> Chain 1 Iteration:  340 / 4000 [  8%]  (Warmup) 
#> Chain 4 Iteration:  250 / 4000 [  6%]  (Warmup) 
#> Chain 1 Iteration:  350 / 4000 [  8%]  (Warmup) 
#> Chain 1 Iteration:  360 / 4000 [  9%]  (Warmup) 
#> Chain 2 Iteration:  280 / 4000 [  7%]  (Warmup) 
#> Chain 3 Iteration:  170 / 4000 [  4%]  (Warmup) 
#> Chain 4 Iteration:  260 / 4000 [  6%]  (Warmup) 
#> Chain 1 Iteration:  370 / 4000 [  9%]  (Warmup) 
#> Chain 2 Iteration:  290 / 4000 [  7%]  (Warmup) 
#> Chain 1 Iteration:  380 / 4000 [  9%]  (Warmup) 
#> Chain 4 Iteration:  270 / 4000 [  6%]  (Warmup) 
#> Chain 1 Iteration:  390 / 4000 [  9%]  (Warmup) 
#> Chain 2 Iteration:  300 / 4000 [  7%]  (Warmup) 
#> Chain 1 Iteration:  400 / 4000 [ 10%]  (Warmup) 
#> Chain 3 Iteration:  180 / 4000 [  4%]  (Warmup) 
#> Chain 4 Iteration:  280 / 4000 [  7%]  (Warmup) 
#> Chain 2 Iteration:  310 / 4000 [  7%]  (Warmup) 
#> Chain 1 Iteration:  410 / 4000 [ 10%]  (Warmup) 
#> Chain 4 Iteration:  290 / 4000 [  7%]  (Warmup) 
#> Chain 3 Iteration:  190 / 4000 [  4%]  (Warmup) 
#> Chain 1 Iteration:  420 / 4000 [ 10%]  (Warmup) 
#> Chain 2 Iteration:  320 / 4000 [  8%]  (Warmup) 
#> Chain 1 Iteration:  430 / 4000 [ 10%]  (Warmup) 
#> Chain 4 Iteration:  300 / 4000 [  7%]  (Warmup) 
#> Chain 2 Iteration:  330 / 4000 [  8%]  (Warmup) 
#> Chain 3 Iteration:  200 / 4000 [  5%]  (Warmup) 
#> Chain 1 Iteration:  440 / 4000 [ 11%]  (Warmup) 
#> Chain 2 Iteration:  340 / 4000 [  8%]  (Warmup) 
#> Chain 3 Iteration:  210 / 4000 [  5%]  (Warmup) 
#> Chain 4 Iteration:  310 / 4000 [  7%]  (Warmup) 
#> Chain 1 Iteration:  450 / 4000 [ 11%]  (Warmup) 
#> Chain 2 Iteration:  350 / 4000 [  8%]  (Warmup) 
#> Chain 4 Iteration:  320 / 4000 [  8%]  (Warmup) 
#> Chain 3 Iteration:  220 / 4000 [  5%]  (Warmup) 
#> Chain 2 Iteration:  360 / 4000 [  9%]  (Warmup) 
#> Chain 4 Iteration:  330 / 4000 [  8%]  (Warmup) 
#> Chain 3 Iteration:  230 / 4000 [  5%]  (Warmup) 
#> Chain 2 Iteration:  370 / 4000 [  9%]  (Warmup) 
#> Chain 3 Iteration:  240 / 4000 [  6%]  (Warmup) 
#> Chain 2 Iteration:  380 / 4000 [  9%]  (Warmup) 
#> Chain 4 Iteration:  340 / 4000 [  8%]  (Warmup) 
#> Chain 3 Iteration:  250 / 4000 [  6%]  (Warmup) 
#> Chain 2 Iteration:  390 / 4000 [  9%]  (Warmup) 
#> Chain 4 Iteration:  350 / 4000 [  8%]  (Warmup) 
#> Chain 2 Iteration:  400 / 4000 [ 10%]  (Warmup) 
#> Chain 2 Iteration:  410 / 4000 [ 10%]  (Warmup) 
#> Chain 1 Iteration:  460 / 4000 [ 11%]  (Warmup) 
#> Chain 4 Iteration:  360 / 4000 [  9%]  (Warmup) 
#> Chain 2 Iteration:  420 / 4000 [ 10%]  (Warmup) 
#> Chain 4 Iteration:  370 / 4000 [  9%]  (Warmup) 
#> Chain 2 Iteration:  430 / 4000 [ 10%]  (Warmup) 
#> Chain 3 Iteration:  260 / 4000 [  6%]  (Warmup) 
#> Chain 2 Iteration:  440 / 4000 [ 11%]  (Warmup) 
#> Chain 4 Iteration:  380 / 4000 [  9%]  (Warmup) 
#> Chain 4 Iteration:  390 / 4000 [  9%]  (Warmup) 
#> Chain 2 Iteration:  450 / 4000 [ 11%]  (Warmup) 
#> Chain 4 Iteration:  400 / 4000 [ 10%]  (Warmup) 
#> Chain 4 Iteration:  410 / 4000 [ 10%]  (Warmup) 
#> Chain 2 Iteration:  460 / 4000 [ 11%]  (Warmup) 
#> Chain 1 Iteration:  470 / 4000 [ 11%]  (Warmup) 
#> Chain 4 Iteration:  420 / 4000 [ 10%]  (Warmup) 
#> Chain 3 Iteration:  270 / 4000 [  6%]  (Warmup) 
#> Chain 4 Iteration:  430 / 4000 [ 10%]  (Warmup) 
#> Chain 4 Iteration:  440 / 4000 [ 11%]  (Warmup) 
#> Chain 4 Iteration:  450 / 4000 [ 11%]  (Warmup) 
#> Chain 2 Iteration:  470 / 4000 [ 11%]  (Warmup) 
#> Chain 3 Iteration:  280 / 4000 [  7%]  (Warmup) 
#> Chain 1 Iteration:  480 / 4000 [ 12%]  (Warmup) 
#> Chain 3 Iteration:  290 / 4000 [  7%]  (Warmup) 
#> Chain 2 Iteration:  480 / 4000 [ 12%]  (Warmup) 
#> Chain 1 Iteration:  490 / 4000 [ 12%]  (Warmup) 
#> Chain 3 Iteration:  300 / 4000 [  7%]  (Warmup) 
#> Chain 4 Iteration:  460 / 4000 [ 11%]  (Warmup) 
#> Chain 2 Iteration:  490 / 4000 [ 12%]  (Warmup) 
#> Chain 3 Iteration:  310 / 4000 [  7%]  (Warmup) 
#> Chain 2 Iteration:  500 / 4000 [ 12%]  (Warmup) 
#> Chain 1 Iteration:  500 / 4000 [ 12%]  (Warmup) 
#> Chain 3 Iteration:  320 / 4000 [  8%]  (Warmup) 
#> Chain 2 Iteration:  510 / 4000 [ 12%]  (Warmup) 
#> Chain 3 Iteration:  330 / 4000 [  8%]  (Warmup) 
#> Chain 1 Iteration:  510 / 4000 [ 12%]  (Warmup) 
#> Chain 2 Iteration:  520 / 4000 [ 13%]  (Warmup) 
#> Chain 3 Iteration:  340 / 4000 [  8%]  (Warmup) 
#> Chain 1 Iteration:  520 / 4000 [ 13%]  (Warmup) 
#> Chain 3 Iteration:  350 / 4000 [  8%]  (Warmup) 
#> Chain 1 Iteration:  530 / 4000 [ 13%]  (Warmup) 
#> Chain 2 Iteration:  530 / 4000 [ 13%]  (Warmup) 
#> Chain 3 Iteration:  360 / 4000 [  9%]  (Warmup) 
#> Chain 1 Iteration:  540 / 4000 [ 13%]  (Warmup) 
#> Chain 2 Iteration:  540 / 4000 [ 13%]  (Warmup) 
#> Chain 3 Iteration:  370 / 4000 [  9%]  (Warmup) 
#> Chain 2 Iteration:  550 / 4000 [ 13%]  (Warmup) 
#> Chain 4 Iteration:  470 / 4000 [ 11%]  (Warmup) 
#> Chain 1 Iteration:  550 / 4000 [ 13%]  (Warmup) 
#> Chain 3 Iteration:  380 / 4000 [  9%]  (Warmup) 
#> Chain 2 Iteration:  560 / 4000 [ 14%]  (Warmup) 
#> Chain 1 Iteration:  560 / 4000 [ 14%]  (Warmup) 
#> Chain 3 Iteration:  390 / 4000 [  9%]  (Warmup) 
#> Chain 2 Iteration:  570 / 4000 [ 14%]  (Warmup) 
#> Chain 1 Iteration:  570 / 4000 [ 14%]  (Warmup) 
#> Chain 3 Iteration:  400 / 4000 [ 10%]  (Warmup) 
#> Chain 2 Iteration:  580 / 4000 [ 14%]  (Warmup) 
#> Chain 1 Iteration:  580 / 4000 [ 14%]  (Warmup) 
#> Chain 3 Iteration:  410 / 4000 [ 10%]  (Warmup) 
#> Chain 2 Iteration:  590 / 4000 [ 14%]  (Warmup) 
#> Chain 4 Iteration:  480 / 4000 [ 12%]  (Warmup) 
#> Chain 1 Iteration:  590 / 4000 [ 14%]  (Warmup) 
#> Chain 3 Iteration:  420 / 4000 [ 10%]  (Warmup) 
#> Chain 2 Iteration:  600 / 4000 [ 15%]  (Warmup) 
#> Chain 3 Iteration:  430 / 4000 [ 10%]  (Warmup) 
#> Chain 1 Iteration:  600 / 4000 [ 15%]  (Warmup) 
#> Chain 2 Iteration:  610 / 4000 [ 15%]  (Warmup) 
#> Chain 4 Iteration:  490 / 4000 [ 12%]  (Warmup) 
#> Chain 2 Iteration:  620 / 4000 [ 15%]  (Warmup) 
#> Chain 3 Iteration:  440 / 4000 [ 11%]  (Warmup) 
#> Chain 1 Iteration:  610 / 4000 [ 15%]  (Warmup) 
#> Chain 2 Iteration:  630 / 4000 [ 15%]  (Warmup) 
#> Chain 3 Iteration:  450 / 4000 [ 11%]  (Warmup) 
#> Chain 4 Iteration:  500 / 4000 [ 12%]  (Warmup) 
#> Chain 1 Iteration:  620 / 4000 [ 15%]  (Warmup) 
#> Chain 2 Iteration:  640 / 4000 [ 16%]  (Warmup) 
#> Chain 4 Iteration:  510 / 4000 [ 12%]  (Warmup) 
#> Chain 1 Iteration:  630 / 4000 [ 15%]  (Warmup) 
#> Chain 2 Iteration:  650 / 4000 [ 16%]  (Warmup) 
#> Chain 2 Iteration:  660 / 4000 [ 16%]  (Warmup) 
#> Chain 1 Iteration:  640 / 4000 [ 16%]  (Warmup) 
#> Chain 4 Iteration:  520 / 4000 [ 13%]  (Warmup) 
#> Chain 2 Iteration:  670 / 4000 [ 16%]  (Warmup) 
#> Chain 1 Iteration:  650 / 4000 [ 16%]  (Warmup) 
#> Chain 2 Iteration:  680 / 4000 [ 17%]  (Warmup) 
#> Chain 4 Iteration:  530 / 4000 [ 13%]  (Warmup) 
#> Chain 2 Iteration:  690 / 4000 [ 17%]  (Warmup) 
#> Chain 1 Iteration:  660 / 4000 [ 16%]  (Warmup) 
#> Chain 2 Iteration:  700 / 4000 [ 17%]  (Warmup) 
#> Chain 4 Iteration:  540 / 4000 [ 13%]  (Warmup) 
#> Chain 3 Iteration:  460 / 4000 [ 11%]  (Warmup) 
#> Chain 1 Iteration:  670 / 4000 [ 16%]  (Warmup) 
#> Chain 2 Iteration:  710 / 4000 [ 17%]  (Warmup) 
#> Chain 4 Iteration:  550 / 4000 [ 13%]  (Warmup) 
#> Chain 2 Iteration:  720 / 4000 [ 18%]  (Warmup) 
#> Chain 4 Iteration:  560 / 4000 [ 14%]  (Warmup) 
#> Chain 1 Iteration:  680 / 4000 [ 17%]  (Warmup) 
#> Chain 2 Iteration:  730 / 4000 [ 18%]  (Warmup) 
#> Chain 4 Iteration:  570 / 4000 [ 14%]  (Warmup) 
#> Chain 1 Iteration:  690 / 4000 [ 17%]  (Warmup) 
#> Chain 2 Iteration:  740 / 4000 [ 18%]  (Warmup) 
#> Chain 4 Iteration:  580 / 4000 [ 14%]  (Warmup) 
#> Chain 2 Iteration:  750 / 4000 [ 18%]  (Warmup) 
#> Chain 1 Iteration:  700 / 4000 [ 17%]  (Warmup) 
#> Chain 2 Iteration:  760 / 4000 [ 19%]  (Warmup) 
#> Chain 4 Iteration:  590 / 4000 [ 14%]  (Warmup) 
#> Chain 4 Iteration:  600 / 4000 [ 15%]  (Warmup) 
#> Chain 2 Iteration:  770 / 4000 [ 19%]  (Warmup) 
#> Chain 1 Iteration:  710 / 4000 [ 17%]  (Warmup) 
#> Chain 4 Iteration:  610 / 4000 [ 15%]  (Warmup) 
#> Chain 2 Iteration:  780 / 4000 [ 19%]  (Warmup) 
#> Chain 1 Iteration:  720 / 4000 [ 18%]  (Warmup) 
#> Chain 4 Iteration:  620 / 4000 [ 15%]  (Warmup) 
#> Chain 2 Iteration:  790 / 4000 [ 19%]  (Warmup) 
#> Chain 1 Iteration:  730 / 4000 [ 18%]  (Warmup) 
#> Chain 4 Iteration:  630 / 4000 [ 15%]  (Warmup) 
#> Chain 2 Iteration:  800 / 4000 [ 20%]  (Warmup) 
#> Chain 1 Iteration:  740 / 4000 [ 18%]  (Warmup) 
#> Chain 4 Iteration:  640 / 4000 [ 16%]  (Warmup) 
#> Chain 2 Iteration:  810 / 4000 [ 20%]  (Warmup) 
#> Chain 1 Iteration:  750 / 4000 [ 18%]  (Warmup) 
#> Chain 3 Iteration:  470 / 4000 [ 11%]  (Warmup) 
#> Chain 4 Iteration:  650 / 4000 [ 16%]  (Warmup) 
#> Chain 2 Iteration:  820 / 4000 [ 20%]  (Warmup) 
#> Chain 1 Iteration:  760 / 4000 [ 19%]  (Warmup) 
#> Chain 2 Iteration:  830 / 4000 [ 20%]  (Warmup) 
#> Chain 4 Iteration:  660 / 4000 [ 16%]  (Warmup) 
#> Chain 1 Iteration:  770 / 4000 [ 19%]  (Warmup) 
#> Chain 2 Iteration:  840 / 4000 [ 21%]  (Warmup) 
#> Chain 4 Iteration:  670 / 4000 [ 16%]  (Warmup) 
#> Chain 1 Iteration:  780 / 4000 [ 19%]  (Warmup) 
#> Chain 2 Iteration:  850 / 4000 [ 21%]  (Warmup) 
#> Chain 4 Iteration:  680 / 4000 [ 17%]  (Warmup) 
#> Chain 1 Iteration:  790 / 4000 [ 19%]  (Warmup) 
#> Chain 4 Iteration:  690 / 4000 [ 17%]  (Warmup) 
#> Chain 1 Iteration:  800 / 4000 [ 20%]  (Warmup) 
#> Chain 3 Iteration:  480 / 4000 [ 12%]  (Warmup) 
#> Chain 4 Iteration:  700 / 4000 [ 17%]  (Warmup) 
#> Chain 1 Iteration:  810 / 4000 [ 20%]  (Warmup) 
#> Chain 2 Iteration:  860 / 4000 [ 21%]  (Warmup) 
#> Chain 4 Iteration:  710 / 4000 [ 17%]  (Warmup) 
#> Chain 1 Iteration:  820 / 4000 [ 20%]  (Warmup) 
#> Chain 4 Iteration:  720 / 4000 [ 18%]  (Warmup) 
#> Chain 3 Iteration:  490 / 4000 [ 12%]  (Warmup) 
#> Chain 4 Iteration:  730 / 4000 [ 18%]  (Warmup) 
#> Chain 1 Iteration:  830 / 4000 [ 20%]  (Warmup) 
#> Chain 3 Iteration:  500 / 4000 [ 12%]  (Warmup) 
#> Chain 4 Iteration:  740 / 4000 [ 18%]  (Warmup) 
#> Chain 1 Iteration:  840 / 4000 [ 21%]  (Warmup) 
#> Chain 1 Iteration:  850 / 4000 [ 21%]  (Warmup) 
#> Chain 2 Iteration:  870 / 4000 [ 21%]  (Warmup) 
#> Chain 3 Iteration:  510 / 4000 [ 12%]  (Warmup) 
#> Chain 4 Iteration:  750 / 4000 [ 18%]  (Warmup) 
#> Chain 3 Iteration:  520 / 4000 [ 13%]  (Warmup) 
#> Chain 3 Iteration:  530 / 4000 [ 13%]  (Warmup) 
#> Chain 4 Iteration:  760 / 4000 [ 19%]  (Warmup) 
#> Chain 2 Iteration:  880 / 4000 [ 22%]  (Warmup) 
#> Chain 4 Iteration:  770 / 4000 [ 19%]  (Warmup) 
#> Chain 3 Iteration:  540 / 4000 [ 13%]  (Warmup) 
#> Chain 1 Iteration:  860 / 4000 [ 21%]  (Warmup) 
#> Chain 2 Iteration:  890 / 4000 [ 22%]  (Warmup) 
#> Chain 4 Iteration:  780 / 4000 [ 19%]  (Warmup) 
#> Chain 3 Iteration:  550 / 4000 [ 13%]  (Warmup) 
#> Chain 2 Iteration:  900 / 4000 [ 22%]  (Warmup) 
#> Chain 4 Iteration:  790 / 4000 [ 19%]  (Warmup) 
#> Chain 3 Iteration:  560 / 4000 [ 14%]  (Warmup) 
#> Chain 2 Iteration:  910 / 4000 [ 22%]  (Warmup) 
#> Chain 4 Iteration:  800 / 4000 [ 20%]  (Warmup) 
#> Chain 2 Iteration:  920 / 4000 [ 23%]  (Warmup) 
#> Chain 3 Iteration:  570 / 4000 [ 14%]  (Warmup) 
#> Chain 4 Iteration:  810 / 4000 [ 20%]  (Warmup) 
#> Chain 2 Iteration:  930 / 4000 [ 23%]  (Warmup) 
#> Chain 1 Iteration:  870 / 4000 [ 21%]  (Warmup) 
#> Chain 4 Iteration:  820 / 4000 [ 20%]  (Warmup) 
#> Chain 3 Iteration:  580 / 4000 [ 14%]  (Warmup) 
#> Chain 2 Iteration:  940 / 4000 [ 23%]  (Warmup) 
#> Chain 4 Iteration:  830 / 4000 [ 20%]  (Warmup) 
#> Chain 2 Iteration:  950 / 4000 [ 23%]  (Warmup) 
#> Chain 3 Iteration:  590 / 4000 [ 14%]  (Warmup) 
#> Chain 4 Iteration:  840 / 4000 [ 21%]  (Warmup) 
#> Chain 2 Iteration:  960 / 4000 [ 24%]  (Warmup) 
#> Chain 1 Iteration:  880 / 4000 [ 22%]  (Warmup) 
#> Chain 4 Iteration:  850 / 4000 [ 21%]  (Warmup) 
#> Chain 3 Iteration:  600 / 4000 [ 15%]  (Warmup) 
#> Chain 2 Iteration:  970 / 4000 [ 24%]  (Warmup) 
#> Chain 2 Iteration:  980 / 4000 [ 24%]  (Warmup) 
#> Chain 1 Iteration:  890 / 4000 [ 22%]  (Warmup) 
#> Chain 3 Iteration:  610 / 4000 [ 15%]  (Warmup) 
#> Chain 2 Iteration:  990 / 4000 [ 24%]  (Warmup) 
#> Chain 1 Iteration:  900 / 4000 [ 22%]  (Warmup) 
#> Chain 3 Iteration:  620 / 4000 [ 15%]  (Warmup) 
#> Chain 2 Iteration: 1000 / 4000 [ 25%]  (Warmup) 
#> Chain 4 Iteration:  860 / 4000 [ 21%]  (Warmup) 
#> Chain 3 Iteration:  630 / 4000 [ 15%]  (Warmup) 
#> Chain 1 Iteration:  910 / 4000 [ 22%]  (Warmup) 
#> Chain 2 Iteration: 1010 / 4000 [ 25%]  (Warmup) 
#> Chain 3 Iteration:  640 / 4000 [ 16%]  (Warmup) 
#> Chain 1 Iteration:  920 / 4000 [ 23%]  (Warmup) 
#> Chain 2 Iteration: 1020 / 4000 [ 25%]  (Warmup) 
#> Chain 3 Iteration:  650 / 4000 [ 16%]  (Warmup) 
#> Chain 2 Iteration: 1030 / 4000 [ 25%]  (Warmup) 
#> Chain 1 Iteration:  930 / 4000 [ 23%]  (Warmup) 
#> Chain 3 Iteration:  660 / 4000 [ 16%]  (Warmup) 
#> Chain 2 Iteration: 1040 / 4000 [ 26%]  (Warmup) 
#> Chain 1 Iteration:  940 / 4000 [ 23%]  (Warmup) 
#> Chain 3 Iteration:  670 / 4000 [ 16%]  (Warmup) 
#> Chain 1 Iteration:  950 / 4000 [ 23%]  (Warmup) 
#> Chain 2 Iteration: 1050 / 4000 [ 26%]  (Warmup) 
#> Chain 3 Iteration:  680 / 4000 [ 17%]  (Warmup) 
#> Chain 4 Iteration:  870 / 4000 [ 21%]  (Warmup) 
#> Chain 1 Iteration:  960 / 4000 [ 24%]  (Warmup) 
#> Chain 3 Iteration:  690 / 4000 [ 17%]  (Warmup) 
#> Chain 2 Iteration: 1060 / 4000 [ 26%]  (Warmup) 
#> Chain 1 Iteration:  970 / 4000 [ 24%]  (Warmup) 
#> Chain 2 Iteration: 1070 / 4000 [ 26%]  (Warmup) 
#> Chain 1 Iteration:  980 / 4000 [ 24%]  (Warmup) 
#> Chain 3 Iteration:  700 / 4000 [ 17%]  (Warmup) 
#> Chain 2 Iteration: 1080 / 4000 [ 27%]  (Warmup) 
#> Chain 1 Iteration:  990 / 4000 [ 24%]  (Warmup) 
#> Chain 4 Iteration:  880 / 4000 [ 22%]  (Warmup) 
#> Chain 2 Iteration: 1090 / 4000 [ 27%]  (Warmup) 
#> Chain 3 Iteration:  710 / 4000 [ 17%]  (Warmup) 
#> Chain 1 Iteration: 1000 / 4000 [ 25%]  (Warmup) 
#> Chain 2 Iteration: 1100 / 4000 [ 27%]  (Warmup) 
#> Chain 4 Iteration:  890 / 4000 [ 22%]  (Warmup) 
#> Chain 1 Iteration: 1010 / 4000 [ 25%]  (Warmup) 
#> Chain 2 Iteration: 1110 / 4000 [ 27%]  (Warmup) 
#> Chain 3 Iteration:  720 / 4000 [ 18%]  (Warmup) 
#> Chain 4 Iteration:  900 / 4000 [ 22%]  (Warmup) 
#> Chain 1 Iteration: 1020 / 4000 [ 25%]  (Warmup) 
#> Chain 2 Iteration: 1120 / 4000 [ 28%]  (Warmup) 
#> Chain 3 Iteration:  730 / 4000 [ 18%]  (Warmup) 
#> Chain 1 Iteration: 1030 / 4000 [ 25%]  (Warmup) 
#> Chain 4 Iteration:  910 / 4000 [ 22%]  (Warmup) 
#> Chain 2 Iteration: 1130 / 4000 [ 28%]  (Warmup) 
#> Chain 1 Iteration: 1040 / 4000 [ 26%]  (Warmup) 
#> Chain 4 Iteration:  920 / 4000 [ 23%]  (Warmup) 
#> Chain 2 Iteration: 1140 / 4000 [ 28%]  (Warmup) 
#> Chain 3 Iteration:  740 / 4000 [ 18%]  (Warmup) 
#> Chain 1 Iteration: 1050 / 4000 [ 26%]  (Warmup) 
#> Chain 4 Iteration:  930 / 4000 [ 23%]  (Warmup) 
#> Chain 2 Iteration: 1150 / 4000 [ 28%]  (Warmup) 
#> Chain 3 Iteration:  750 / 4000 [ 18%]  (Warmup) 
#> Chain 1 Iteration: 1060 / 4000 [ 26%]  (Warmup) 
#> Chain 4 Iteration:  940 / 4000 [ 23%]  (Warmup) 
#> Chain 2 Iteration: 1160 / 4000 [ 29%]  (Warmup) 
#> Chain 3 Iteration:  760 / 4000 [ 19%]  (Warmup) 
#> Chain 1 Iteration: 1070 / 4000 [ 26%]  (Warmup) 
#> Chain 4 Iteration:  950 / 4000 [ 23%]  (Warmup) 
#> Chain 2 Iteration: 1170 / 4000 [ 29%]  (Warmup) 
#> Chain 1 Iteration: 1080 / 4000 [ 27%]  (Warmup) 
#> Chain 4 Iteration:  960 / 4000 [ 24%]  (Warmup) 
#> Chain 3 Iteration:  770 / 4000 [ 19%]  (Warmup) 
#> Chain 2 Iteration: 1180 / 4000 [ 29%]  (Warmup) 
#> Chain 4 Iteration:  970 / 4000 [ 24%]  (Warmup) 
#> Chain 1 Iteration: 1090 / 4000 [ 27%]  (Warmup) 
#> Chain 2 Iteration: 1190 / 4000 [ 29%]  (Warmup) 
#> Chain 3 Iteration:  780 / 4000 [ 19%]  (Warmup) 
#> Chain 4 Iteration:  980 / 4000 [ 24%]  (Warmup) 
#> Chain 2 Iteration: 1200 / 4000 [ 30%]  (Warmup) 
#> Chain 1 Iteration: 1100 / 4000 [ 27%]  (Warmup) 
#> Chain 4 Iteration:  990 / 4000 [ 24%]  (Warmup) 
#> Chain 2 Iteration: 1210 / 4000 [ 30%]  (Warmup) 
#> Chain 1 Iteration: 1110 / 4000 [ 27%]  (Warmup) 
#> Chain 3 Iteration:  790 / 4000 [ 19%]  (Warmup) 
#> Chain 4 Iteration: 1000 / 4000 [ 25%]  (Warmup) 
#> Chain 2 Iteration: 1220 / 4000 [ 30%]  (Warmup) 
#> Chain 4 Iteration: 1010 / 4000 [ 25%]  (Warmup) 
#> Chain 1 Iteration: 1120 / 4000 [ 28%]  (Warmup) 
#> Chain 3 Iteration:  800 / 4000 [ 20%]  (Warmup) 
#> Chain 4 Iteration: 1020 / 4000 [ 25%]  (Warmup) 
#> Chain 1 Iteration: 1130 / 4000 [ 28%]  (Warmup) 
#> Chain 2 Iteration: 1230 / 4000 [ 30%]  (Warmup) 
#> Chain 4 Iteration: 1030 / 4000 [ 25%]  (Warmup) 
#> Chain 3 Iteration:  810 / 4000 [ 20%]  (Warmup) 
#> Chain 2 Iteration: 1240 / 4000 [ 31%]  (Warmup) 
#> Chain 1 Iteration: 1140 / 4000 [ 28%]  (Warmup) 
#> Chain 4 Iteration: 1040 / 4000 [ 26%]  (Warmup) 
#> Chain 2 Iteration: 1250 / 4000 [ 31%]  (Warmup) 
#> Chain 3 Iteration:  820 / 4000 [ 20%]  (Warmup) 
#> Chain 1 Iteration: 1150 / 4000 [ 28%]  (Warmup) 
#> Chain 2 Iteration: 1260 / 4000 [ 31%]  (Warmup) 
#> Chain 4 Iteration: 1050 / 4000 [ 26%]  (Warmup) 
#> Chain 2 Iteration: 1270 / 4000 [ 31%]  (Warmup) 
#> Chain 3 Iteration:  830 / 4000 [ 20%]  (Warmup) 
#> Chain 1 Iteration: 1160 / 4000 [ 29%]  (Warmup) 
#> Chain 4 Iteration: 1060 / 4000 [ 26%]  (Warmup) 
#> Chain 2 Iteration: 1280 / 4000 [ 32%]  (Warmup) 
#> Chain 1 Iteration: 1170 / 4000 [ 29%]  (Warmup) 
#> Chain 4 Iteration: 1070 / 4000 [ 26%]  (Warmup) 
#> Chain 3 Iteration:  840 / 4000 [ 21%]  (Warmup) 
#> Chain 2 Iteration: 1290 / 4000 [ 32%]  (Warmup) 
#> Chain 1 Iteration: 1180 / 4000 [ 29%]  (Warmup) 
#> Chain 2 Iteration: 1300 / 4000 [ 32%]  (Warmup) 
#> Chain 3 Iteration:  850 / 4000 [ 21%]  (Warmup) 
#> Chain 1 Iteration: 1190 / 4000 [ 29%]  (Warmup) 
#> Chain 4 Iteration: 1080 / 4000 [ 27%]  (Warmup) 
#> Chain 2 Iteration: 1310 / 4000 [ 32%]  (Warmup) 
#> Chain 3 Iteration:  860 / 4000 [ 21%]  (Warmup) 
#> Chain 1 Iteration: 1200 / 4000 [ 30%]  (Warmup) 
#> Chain 4 Iteration: 1090 / 4000 [ 27%]  (Warmup) 
#> Chain 1 Iteration: 1210 / 4000 [ 30%]  (Warmup) 
#> Chain 2 Iteration: 1320 / 4000 [ 33%]  (Warmup) 
#> Chain 3 Iteration:  870 / 4000 [ 21%]  (Warmup) 
#> Chain 4 Iteration: 1100 / 4000 [ 27%]  (Warmup) 
#> Chain 1 Iteration: 1220 / 4000 [ 30%]  (Warmup) 
#> Chain 2 Iteration: 1330 / 4000 [ 33%]  (Warmup) 
#> Chain 4 Iteration: 1110 / 4000 [ 27%]  (Warmup) 
#> Chain 3 Iteration:  880 / 4000 [ 22%]  (Warmup) 
#> Chain 1 Iteration: 1230 / 4000 [ 30%]  (Warmup) 
#> Chain 2 Iteration: 1340 / 4000 [ 33%]  (Warmup) 
#> Chain 4 Iteration: 1120 / 4000 [ 28%]  (Warmup) 
#> Chain 3 Iteration:  890 / 4000 [ 22%]  (Warmup) 
#> Chain 1 Iteration: 1240 / 4000 [ 31%]  (Warmup) 
#> Chain 2 Iteration: 1350 / 4000 [ 33%]  (Warmup) 
#> Chain 4 Iteration: 1130 / 4000 [ 28%]  (Warmup) 
#> Chain 3 Iteration:  900 / 4000 [ 22%]  (Warmup) 
#> Chain 1 Iteration: 1250 / 4000 [ 31%]  (Warmup) 
#> Chain 2 Iteration: 1360 / 4000 [ 34%]  (Warmup) 
#> Chain 4 Iteration: 1140 / 4000 [ 28%]  (Warmup) 
#> Chain 3 Iteration:  910 / 4000 [ 22%]  (Warmup) 
#> Chain 1 Iteration: 1260 / 4000 [ 31%]  (Warmup) 
#> Chain 2 Iteration: 1370 / 4000 [ 34%]  (Warmup) 
#> Chain 4 Iteration: 1150 / 4000 [ 28%]  (Warmup) 
#> Chain 3 Iteration:  920 / 4000 [ 23%]  (Warmup) 
#> Chain 2 Iteration: 1380 / 4000 [ 34%]  (Warmup) 
#> Chain 4 Iteration: 1160 / 4000 [ 29%]  (Warmup) 
#> Chain 1 Iteration: 1270 / 4000 [ 31%]  (Warmup) 
#> Chain 3 Iteration:  930 / 4000 [ 23%]  (Warmup) 
#> Chain 2 Iteration: 1390 / 4000 [ 34%]  (Warmup) 
#> Chain 4 Iteration: 1170 / 4000 [ 29%]  (Warmup) 
#> Chain 2 Iteration: 1400 / 4000 [ 35%]  (Warmup) 
#> Chain 3 Iteration:  940 / 4000 [ 23%]  (Warmup) 
#> Chain 1 Iteration: 1280 / 4000 [ 32%]  (Warmup) 
#> Chain 4 Iteration: 1180 / 4000 [ 29%]  (Warmup) 
#> Chain 2 Iteration: 1410 / 4000 [ 35%]  (Warmup) 
#> Chain 3 Iteration:  950 / 4000 [ 23%]  (Warmup) 
#> Chain 2 Iteration: 1420 / 4000 [ 35%]  (Warmup) 
#> Chain 4 Iteration: 1190 / 4000 [ 29%]  (Warmup) 
#> Chain 1 Iteration: 1290 / 4000 [ 32%]  (Warmup) 
#> Chain 3 Iteration:  960 / 4000 [ 24%]  (Warmup) 
#> Chain 4 Iteration: 1200 / 4000 [ 30%]  (Warmup) 
#> Chain 3 Iteration:  970 / 4000 [ 24%]  (Warmup) 
#> Chain 2 Iteration: 1430 / 4000 [ 35%]  (Warmup) 
#> Chain 1 Iteration: 1300 / 4000 [ 32%]  (Warmup) 
#> Chain 4 Iteration: 1210 / 4000 [ 30%]  (Warmup) 
#> Chain 3 Iteration:  980 / 4000 [ 24%]  (Warmup) 
#> Chain 1 Iteration: 1310 / 4000 [ 32%]  (Warmup) 
#> Chain 2 Iteration: 1440 / 4000 [ 36%]  (Warmup) 
#> Chain 3 Iteration:  990 / 4000 [ 24%]  (Warmup) 
#> Chain 4 Iteration: 1220 / 4000 [ 30%]  (Warmup) 
#> Chain 3 Iteration: 1000 / 4000 [ 25%]  (Warmup) 
#> Chain 2 Iteration: 1450 / 4000 [ 36%]  (Warmup) 
#> Chain 1 Iteration: 1320 / 4000 [ 33%]  (Warmup) 
#> Chain 4 Iteration: 1230 / 4000 [ 30%]  (Warmup) 
#> Chain 3 Iteration: 1010 / 4000 [ 25%]  (Warmup) 
#> Chain 1 Iteration: 1330 / 4000 [ 33%]  (Warmup) 
#> Chain 4 Iteration: 1240 / 4000 [ 31%]  (Warmup) 
#> Chain 2 Iteration: 1460 / 4000 [ 36%]  (Warmup) 
#> Chain 3 Iteration: 1020 / 4000 [ 25%]  (Warmup) 
#> Chain 1 Iteration: 1340 / 4000 [ 33%]  (Warmup) 
#> Chain 3 Iteration: 1030 / 4000 [ 25%]  (Warmup) 
#> Chain 4 Iteration: 1250 / 4000 [ 31%]  (Warmup) 
#> Chain 2 Iteration: 1470 / 4000 [ 36%]  (Warmup) 
#> Chain 3 Iteration: 1040 / 4000 [ 26%]  (Warmup) 
#> Chain 1 Iteration: 1350 / 4000 [ 33%]  (Warmup) 
#> Chain 2 Iteration: 1480 / 4000 [ 37%]  (Warmup) 
#> Chain 4 Iteration: 1260 / 4000 [ 31%]  (Warmup) 
#> Chain 3 Iteration: 1050 / 4000 [ 26%]  (Warmup) 
#> Chain 1 Iteration: 1360 / 4000 [ 34%]  (Warmup) 
#> Chain 3 Iteration: 1060 / 4000 [ 26%]  (Warmup) 
#> Chain 2 Iteration: 1490 / 4000 [ 37%]  (Warmup) 
#> Chain 4 Iteration: 1270 / 4000 [ 31%]  (Warmup) 
#> Chain 1 Iteration: 1370 / 4000 [ 34%]  (Warmup) 
#> Chain 3 Iteration: 1070 / 4000 [ 26%]  (Warmup) 
#> Chain 2 Iteration: 1500 / 4000 [ 37%]  (Warmup) 
#> Chain 3 Iteration: 1080 / 4000 [ 27%]  (Warmup) 
#> Chain 1 Iteration: 1380 / 4000 [ 34%]  (Warmup) 
#> Chain 4 Iteration: 1280 / 4000 [ 32%]  (Warmup) 
#> Chain 2 Iteration: 1510 / 4000 [ 37%]  (Warmup) 
#> Chain 3 Iteration: 1090 / 4000 [ 27%]  (Warmup) 
#> Chain 4 Iteration: 1290 / 4000 [ 32%]  (Warmup) 
#> Chain 1 Iteration: 1390 / 4000 [ 34%]  (Warmup) 
#> Chain 2 Iteration: 1520 / 4000 [ 38%]  (Warmup) 
#> Chain 3 Iteration: 1100 / 4000 [ 27%]  (Warmup) 
#> Chain 4 Iteration: 1300 / 4000 [ 32%]  (Warmup) 
#> Chain 1 Iteration: 1400 / 4000 [ 35%]  (Warmup) 
#> Chain 2 Iteration: 1530 / 4000 [ 38%]  (Warmup) 
#> Chain 3 Iteration: 1110 / 4000 [ 27%]  (Warmup) 
#> Chain 4 Iteration: 1310 / 4000 [ 32%]  (Warmup) 
#> Chain 1 Iteration: 1410 / 4000 [ 35%]  (Warmup) 
#> Chain 4 Iteration: 1320 / 4000 [ 33%]  (Warmup) 
#> Chain 2 Iteration: 1540 / 4000 [ 38%]  (Warmup) 
#> Chain 3 Iteration: 1120 / 4000 [ 28%]  (Warmup) 
#> Chain 4 Iteration: 1330 / 4000 [ 33%]  (Warmup) 
#> Chain 1 Iteration: 1420 / 4000 [ 35%]  (Warmup) 
#> Chain 2 Iteration: 1550 / 4000 [ 38%]  (Warmup) 
#> Chain 3 Iteration: 1130 / 4000 [ 28%]  (Warmup) 
#> Chain 4 Iteration: 1340 / 4000 [ 33%]  (Warmup) 
#> Chain 2 Iteration: 1560 / 4000 [ 39%]  (Warmup) 
#> Chain 1 Iteration: 1430 / 4000 [ 35%]  (Warmup) 
#> Chain 4 Iteration: 1350 / 4000 [ 33%]  (Warmup) 
#> Chain 3 Iteration: 1140 / 4000 [ 28%]  (Warmup) 
#> Chain 2 Iteration: 1570 / 4000 [ 39%]  (Warmup) 
#> Chain 4 Iteration: 1360 / 4000 [ 34%]  (Warmup) 
#> Chain 1 Iteration: 1440 / 4000 [ 36%]  (Warmup) 
#> Chain 2 Iteration: 1580 / 4000 [ 39%]  (Warmup) 
#> Chain 4 Iteration: 1370 / 4000 [ 34%]  (Warmup) 
#> Chain 3 Iteration: 1150 / 4000 [ 28%]  (Warmup) 
#> Chain 2 Iteration: 1590 / 4000 [ 39%]  (Warmup) 
#> Chain 4 Iteration: 1380 / 4000 [ 34%]  (Warmup) 
#> Chain 1 Iteration: 1450 / 4000 [ 36%]  (Warmup) 
#> Chain 4 Iteration: 1390 / 4000 [ 34%]  (Warmup) 
#> Chain 2 Iteration: 1600 / 4000 [ 40%]  (Warmup) 
#> Chain 3 Iteration: 1160 / 4000 [ 29%]  (Warmup) 
#> Chain 1 Iteration: 1460 / 4000 [ 36%]  (Warmup) 
#> Chain 4 Iteration: 1400 / 4000 [ 35%]  (Warmup) 
#> Chain 2 Iteration: 1610 / 4000 [ 40%]  (Warmup) 
#> Chain 3 Iteration: 1170 / 4000 [ 29%]  (Warmup) 
#> Chain 4 Iteration: 1410 / 4000 [ 35%]  (Warmup) 
#> Chain 2 Iteration: 1620 / 4000 [ 40%]  (Warmup) 
#> Chain 1 Iteration: 1470 / 4000 [ 36%]  (Warmup) 
#> Chain 4 Iteration: 1420 / 4000 [ 35%]  (Warmup) 
#> Chain 2 Iteration: 1630 / 4000 [ 40%]  (Warmup) 
#> Chain 3 Iteration: 1180 / 4000 [ 29%]  (Warmup) 
#> Chain 1 Iteration: 1480 / 4000 [ 37%]  (Warmup) 
#> Chain 4 Iteration: 1430 / 4000 [ 35%]  (Warmup) 
#> Chain 3 Iteration: 1190 / 4000 [ 29%]  (Warmup) 
#> Chain 1 Iteration: 1490 / 4000 [ 37%]  (Warmup) 
#> Chain 2 Iteration: 1640 / 4000 [ 41%]  (Warmup) 
#> Chain 3 Iteration: 1200 / 4000 [ 30%]  (Warmup) 
#> Chain 4 Iteration: 1440 / 4000 [ 36%]  (Warmup) 
#> Chain 2 Iteration: 1650 / 4000 [ 41%]  (Warmup) 
#> Chain 3 Iteration: 1210 / 4000 [ 30%]  (Warmup) 
#> Chain 1 Iteration: 1500 / 4000 [ 37%]  (Warmup) 
#> Chain 4 Iteration: 1450 / 4000 [ 36%]  (Warmup) 
#> Chain 3 Iteration: 1220 / 4000 [ 30%]  (Warmup) 
#> Chain 1 Iteration: 1510 / 4000 [ 37%]  (Warmup) 
#> Chain 2 Iteration: 1660 / 4000 [ 41%]  (Warmup) 
#> Chain 3 Iteration: 1230 / 4000 [ 30%]  (Warmup) 
#> Chain 4 Iteration: 1460 / 4000 [ 36%]  (Warmup) 
#> Chain 1 Iteration: 1520 / 4000 [ 38%]  (Warmup) 
#> Chain 3 Iteration: 1240 / 4000 [ 31%]  (Warmup) 
#> Chain 2 Iteration: 1670 / 4000 [ 41%]  (Warmup) 
#> Chain 4 Iteration: 1470 / 4000 [ 36%]  (Warmup) 
#> Chain 3 Iteration: 1250 / 4000 [ 31%]  (Warmup) 
#> Chain 1 Iteration: 1530 / 4000 [ 38%]  (Warmup) 
#> Chain 2 Iteration: 1680 / 4000 [ 42%]  (Warmup) 
#> Chain 4 Iteration: 1480 / 4000 [ 37%]  (Warmup) 
#> Chain 3 Iteration: 1260 / 4000 [ 31%]  (Warmup) 
#> Chain 1 Iteration: 1540 / 4000 [ 38%]  (Warmup) 
#> Chain 2 Iteration: 1690 / 4000 [ 42%]  (Warmup) 
#> Chain 3 Iteration: 1270 / 4000 [ 31%]  (Warmup) 
#> Chain 4 Iteration: 1490 / 4000 [ 37%]  (Warmup) 
#> Chain 3 Iteration: 1280 / 4000 [ 32%]  (Warmup) 
#> Chain 1 Iteration: 1550 / 4000 [ 38%]  (Warmup) 
#> Chain 2 Iteration: 1700 / 4000 [ 42%]  (Warmup) 
#> Chain 3 Iteration: 1290 / 4000 [ 32%]  (Warmup) 
#> Chain 4 Iteration: 1500 / 4000 [ 37%]  (Warmup) 
#> Chain 1 Iteration: 1560 / 4000 [ 39%]  (Warmup) 
#> Chain 3 Iteration: 1300 / 4000 [ 32%]  (Warmup) 
#> Chain 2 Iteration: 1710 / 4000 [ 42%]  (Warmup) 
#> Chain 3 Iteration: 1310 / 4000 [ 32%]  (Warmup) 
#> Chain 4 Iteration: 1510 / 4000 [ 37%]  (Warmup) 
#> Chain 1 Iteration: 1570 / 4000 [ 39%]  (Warmup) 
#> Chain 3 Iteration: 1320 / 4000 [ 33%]  (Warmup) 
#> Chain 2 Iteration: 1720 / 4000 [ 43%]  (Warmup) 
#> Chain 3 Iteration: 1330 / 4000 [ 33%]  (Warmup) 
#> Chain 4 Iteration: 1520 / 4000 [ 38%]  (Warmup) 
#> Chain 1 Iteration: 1580 / 4000 [ 39%]  (Warmup) 
#> Chain 2 Iteration: 1730 / 4000 [ 43%]  (Warmup) 
#> Chain 3 Iteration: 1340 / 4000 [ 33%]  (Warmup) 
#> Chain 2 Iteration: 1740 / 4000 [ 43%]  (Warmup) 
#> Chain 4 Iteration: 1530 / 4000 [ 38%]  (Warmup) 
#> Chain 1 Iteration: 1590 / 4000 [ 39%]  (Warmup) 
#> Chain 3 Iteration: 1350 / 4000 [ 33%]  (Warmup) 
#> Chain 2 Iteration: 1750 / 4000 [ 43%]  (Warmup) 
#> Chain 4 Iteration: 1540 / 4000 [ 38%]  (Warmup) 
#> Chain 1 Iteration: 1600 / 4000 [ 40%]  (Warmup) 
#> Chain 2 Iteration: 1760 / 4000 [ 44%]  (Warmup) 
#> Chain 3 Iteration: 1360 / 4000 [ 34%]  (Warmup) 
#> Chain 4 Iteration: 1550 / 4000 [ 38%]  (Warmup) 
#> Chain 1 Iteration: 1610 / 4000 [ 40%]  (Warmup) 
#> Chain 2 Iteration: 1770 / 4000 [ 44%]  (Warmup) 
#> Chain 3 Iteration: 1370 / 4000 [ 34%]  (Warmup) 
#> Chain 4 Iteration: 1560 / 4000 [ 39%]  (Warmup) 
#> Chain 3 Iteration: 1380 / 4000 [ 34%]  (Warmup) 
#> Chain 2 Iteration: 1780 / 4000 [ 44%]  (Warmup) 
#> Chain 4 Iteration: 1570 / 4000 [ 39%]  (Warmup) 
#> Chain 1 Iteration: 1620 / 4000 [ 40%]  (Warmup) 
#> Chain 3 Iteration: 1390 / 4000 [ 34%]  (Warmup) 
#> Chain 4 Iteration: 1580 / 4000 [ 39%]  (Warmup) 
#> Chain 2 Iteration: 1790 / 4000 [ 44%]  (Warmup) 
#> Chain 3 Iteration: 1400 / 4000 [ 35%]  (Warmup) 
#> Chain 1 Iteration: 1630 / 4000 [ 40%]  (Warmup) 
#> Chain 3 Iteration: 1410 / 4000 [ 35%]  (Warmup) 
#> Chain 4 Iteration: 1590 / 4000 [ 39%]  (Warmup) 
#> Chain 2 Iteration: 1800 / 4000 [ 45%]  (Warmup) 
#> Chain 4 Iteration: 1600 / 4000 [ 40%]  (Warmup) 
#> Chain 1 Iteration: 1640 / 4000 [ 41%]  (Warmup) 
#> Chain 3 Iteration: 1420 / 4000 [ 35%]  (Warmup) 
#> Chain 2 Iteration: 1810 / 4000 [ 45%]  (Warmup) 
#> Chain 4 Iteration: 1610 / 4000 [ 40%]  (Warmup) 
#> Chain 1 Iteration: 1650 / 4000 [ 41%]  (Warmup) 
#> Chain 3 Iteration: 1430 / 4000 [ 35%]  (Warmup) 
#> Chain 2 Iteration: 1820 / 4000 [ 45%]  (Warmup) 
#> Chain 4 Iteration: 1620 / 4000 [ 40%]  (Warmup) 
#> Chain 3 Iteration: 1440 / 4000 [ 36%]  (Warmup) 
#> Chain 4 Iteration: 1630 / 4000 [ 40%]  (Warmup) 
#> Chain 2 Iteration: 1830 / 4000 [ 45%]  (Warmup) 
#> Chain 1 Iteration: 1660 / 4000 [ 41%]  (Warmup) 
#> Chain 3 Iteration: 1450 / 4000 [ 36%]  (Warmup) 
#> Chain 4 Iteration: 1640 / 4000 [ 41%]  (Warmup) 
#> Chain 2 Iteration: 1840 / 4000 [ 46%]  (Warmup) 
#> Chain 3 Iteration: 1460 / 4000 [ 36%]  (Warmup) 
#> Chain 1 Iteration: 1670 / 4000 [ 41%]  (Warmup) 
#> Chain 4 Iteration: 1650 / 4000 [ 41%]  (Warmup) 
#> Chain 2 Iteration: 1850 / 4000 [ 46%]  (Warmup) 
#> Chain 4 Iteration: 1660 / 4000 [ 41%]  (Warmup) 
#> Chain 3 Iteration: 1470 / 4000 [ 36%]  (Warmup) 
#> Chain 1 Iteration: 1680 / 4000 [ 42%]  (Warmup) 
#> Chain 4 Iteration: 1670 / 4000 [ 41%]  (Warmup) 
#> Chain 3 Iteration: 1480 / 4000 [ 37%]  (Warmup) 
#> Chain 2 Iteration: 1860 / 4000 [ 46%]  (Warmup) 
#> Chain 3 Iteration: 1490 / 4000 [ 37%]  (Warmup) 
#> Chain 4 Iteration: 1680 / 4000 [ 42%]  (Warmup) 
#> Chain 1 Iteration: 1690 / 4000 [ 42%]  (Warmup) 
#> Chain 2 Iteration: 1870 / 4000 [ 46%]  (Warmup) 
#> Chain 3 Iteration: 1500 / 4000 [ 37%]  (Warmup) 
#> Chain 4 Iteration: 1690 / 4000 [ 42%]  (Warmup) 
#> Chain 2 Iteration: 1880 / 4000 [ 47%]  (Warmup) 
#> Chain 3 Iteration: 1510 / 4000 [ 37%]  (Warmup) 
#> Chain 4 Iteration: 1700 / 4000 [ 42%]  (Warmup) 
#> Chain 1 Iteration: 1700 / 4000 [ 42%]  (Warmup) 
#> Chain 2 Iteration: 1890 / 4000 [ 47%]  (Warmup) 
#> Chain 3 Iteration: 1520 / 4000 [ 38%]  (Warmup) 
#> Chain 4 Iteration: 1710 / 4000 [ 42%]  (Warmup) 
#> Chain 1 Iteration: 1710 / 4000 [ 42%]  (Warmup) 
#> Chain 2 Iteration: 1900 / 4000 [ 47%]  (Warmup) 
#> Chain 3 Iteration: 1530 / 4000 [ 38%]  (Warmup) 
#> Chain 2 Iteration: 1910 / 4000 [ 47%]  (Warmup) 
#> Chain 4 Iteration: 1720 / 4000 [ 43%]  (Warmup) 
#> Chain 1 Iteration: 1720 / 4000 [ 43%]  (Warmup) 
#> Chain 3 Iteration: 1540 / 4000 [ 38%]  (Warmup) 
#> Chain 2 Iteration: 1920 / 4000 [ 48%]  (Warmup) 
#> Chain 4 Iteration: 1730 / 4000 [ 43%]  (Warmup) 
#> Chain 1 Iteration: 1730 / 4000 [ 43%]  (Warmup) 
#> Chain 2 Iteration: 1930 / 4000 [ 48%]  (Warmup) 
#> Chain 3 Iteration: 1550 / 4000 [ 38%]  (Warmup) 
#> Chain 4 Iteration: 1740 / 4000 [ 43%]  (Warmup) 
#> Chain 1 Iteration: 1740 / 4000 [ 43%]  (Warmup) 
#> Chain 2 Iteration: 1940 / 4000 [ 48%]  (Warmup) 
#> Chain 3 Iteration: 1560 / 4000 [ 39%]  (Warmup) 
#> Chain 4 Iteration: 1750 / 4000 [ 43%]  (Warmup) 
#> Chain 3 Iteration: 1570 / 4000 [ 39%]  (Warmup) 
#> Chain 1 Iteration: 1750 / 4000 [ 43%]  (Warmup) 
#> Chain 2 Iteration: 1950 / 4000 [ 48%]  (Warmup) 
#> Chain 3 Iteration: 1580 / 4000 [ 39%]  (Warmup) 
#> Chain 4 Iteration: 1760 / 4000 [ 44%]  (Warmup) 
#> Chain 1 Iteration: 1760 / 4000 [ 44%]  (Warmup) 
#> Chain 3 Iteration: 1590 / 4000 [ 39%]  (Warmup) 
#> Chain 4 Iteration: 1770 / 4000 [ 44%]  (Warmup) 
#> Chain 1 Iteration: 1770 / 4000 [ 44%]  (Warmup) 
#> Chain 3 Iteration: 1600 / 4000 [ 40%]  (Warmup) 
#> Chain 2 Iteration: 1960 / 4000 [ 49%]  (Warmup) 
#> Chain 4 Iteration: 1780 / 4000 [ 44%]  (Warmup) 
#> Chain 3 Iteration: 1610 / 4000 [ 40%]  (Warmup) 
#> Chain 1 Iteration: 1780 / 4000 [ 44%]  (Warmup) 
#> Chain 4 Iteration: 1790 / 4000 [ 44%]  (Warmup) 
#> Chain 3 Iteration: 1620 / 4000 [ 40%]  (Warmup) 
#> Chain 1 Iteration: 1790 / 4000 [ 44%]  (Warmup) 
#> Chain 3 Iteration: 1630 / 4000 [ 40%]  (Warmup) 
#> Chain 4 Iteration: 1800 / 4000 [ 45%]  (Warmup) 
#> Chain 1 Iteration: 1800 / 4000 [ 45%]  (Warmup) 
#> Chain 3 Iteration: 1640 / 4000 [ 41%]  (Warmup) 
#> Chain 4 Iteration: 1810 / 4000 [ 45%]  (Warmup) 
#> Chain 2 Iteration: 1970 / 4000 [ 49%]  (Warmup) 
#> Chain 1 Iteration: 1810 / 4000 [ 45%]  (Warmup) 
#> Chain 3 Iteration: 1650 / 4000 [ 41%]  (Warmup) 
#> Chain 4 Iteration: 1820 / 4000 [ 45%]  (Warmup) 
#> Chain 1 Iteration: 1820 / 4000 [ 45%]  (Warmup) 
#> Chain 3 Iteration: 1660 / 4000 [ 41%]  (Warmup) 
#> Chain 3 Iteration: 1670 / 4000 [ 41%]  (Warmup) 
#> Chain 1 Iteration: 1830 / 4000 [ 45%]  (Warmup) 
#> Chain 4 Iteration: 1830 / 4000 [ 45%]  (Warmup) 
#> Chain 2 Iteration: 1980 / 4000 [ 49%]  (Warmup) 
#> Chain 3 Iteration: 1680 / 4000 [ 42%]  (Warmup) 
#> Chain 3 Iteration: 1690 / 4000 [ 42%]  (Warmup) 
#> Chain 4 Iteration: 1840 / 4000 [ 46%]  (Warmup) 
#> Chain 1 Iteration: 1840 / 4000 [ 46%]  (Warmup) 
#> Chain 2 Iteration: 1990 / 4000 [ 49%]  (Warmup) 
#> Chain 3 Iteration: 1700 / 4000 [ 42%]  (Warmup) 
#> Chain 4 Iteration: 1850 / 4000 [ 46%]  (Warmup) 
#> Chain 2 Iteration: 2000 / 4000 [ 50%]  (Warmup) 
#> Chain 2 Iteration: 2001 / 4000 [ 50%]  (Sampling) 
#> Chain 3 Iteration: 1710 / 4000 [ 42%]  (Warmup) 
#> Chain 1 Iteration: 1850 / 4000 [ 46%]  (Warmup) 
#> Chain 3 Iteration: 1720 / 4000 [ 43%]  (Warmup) 
#> Chain 4 Iteration: 1860 / 4000 [ 46%]  (Warmup) 
#> Chain 2 Iteration: 2010 / 4000 [ 50%]  (Sampling) 
#> Chain 1 Iteration: 1860 / 4000 [ 46%]  (Warmup) 
#> Chain 3 Iteration: 1730 / 4000 [ 43%]  (Warmup) 
#> Chain 4 Iteration: 1870 / 4000 [ 46%]  (Warmup) 
#> Chain 3 Iteration: 1740 / 4000 [ 43%]  (Warmup) 
#> Chain 2 Iteration: 2020 / 4000 [ 50%]  (Sampling) 
#> Chain 1 Iteration: 1870 / 4000 [ 46%]  (Warmup) 
#> Chain 4 Iteration: 1880 / 4000 [ 47%]  (Warmup) 
#> Chain 3 Iteration: 1750 / 4000 [ 43%]  (Warmup) 
#> Chain 1 Iteration: 1880 / 4000 [ 47%]  (Warmup) 
#> Chain 2 Iteration: 2030 / 4000 [ 50%]  (Sampling) 
#> Chain 3 Iteration: 1760 / 4000 [ 44%]  (Warmup) 
#> Chain 4 Iteration: 1890 / 4000 [ 47%]  (Warmup) 
#> Chain 3 Iteration: 1770 / 4000 [ 44%]  (Warmup) 
#> Chain 4 Iteration: 1900 / 4000 [ 47%]  (Warmup) 
#> Chain 1 Iteration: 1890 / 4000 [ 47%]  (Warmup) 
#> Chain 2 Iteration: 2040 / 4000 [ 51%]  (Sampling) 
#> Chain 3 Iteration: 1780 / 4000 [ 44%]  (Warmup) 
#> Chain 4 Iteration: 1910 / 4000 [ 47%]  (Warmup) 
#> Chain 1 Iteration: 1900 / 4000 [ 47%]  (Warmup) 
#> Chain 3 Iteration: 1790 / 4000 [ 44%]  (Warmup) 
#> Chain 2 Iteration: 2050 / 4000 [ 51%]  (Sampling) 
#> Chain 4 Iteration: 1920 / 4000 [ 48%]  (Warmup) 
#> Chain 3 Iteration: 1800 / 4000 [ 45%]  (Warmup) 
#> Chain 4 Iteration: 1930 / 4000 [ 48%]  (Warmup) 
#> Chain 1 Iteration: 1910 / 4000 [ 47%]  (Warmup) 
#> Chain 2 Iteration: 2060 / 4000 [ 51%]  (Sampling) 
#> Chain 3 Iteration: 1810 / 4000 [ 45%]  (Warmup) 
#> Chain 4 Iteration: 1940 / 4000 [ 48%]  (Warmup) 
#> Chain 3 Iteration: 1820 / 4000 [ 45%]  (Warmup) 
#> Chain 1 Iteration: 1920 / 4000 [ 48%]  (Warmup) 
#> Chain 4 Iteration: 1950 / 4000 [ 48%]  (Warmup) 
#> Chain 2 Iteration: 2070 / 4000 [ 51%]  (Sampling) 
#> Chain 3 Iteration: 1830 / 4000 [ 45%]  (Warmup) 
#> Chain 3 Iteration: 1840 / 4000 [ 46%]  (Warmup) 
#> Chain 1 Iteration: 1930 / 4000 [ 48%]  (Warmup) 
#> Chain 2 Iteration: 2080 / 4000 [ 52%]  (Sampling) 
#> Chain 3 Iteration: 1850 / 4000 [ 46%]  (Warmup) 
#> Chain 1 Iteration: 1940 / 4000 [ 48%]  (Warmup) 
#> Chain 3 Iteration: 1860 / 4000 [ 46%]  (Warmup) 
#> Chain 2 Iteration: 2090 / 4000 [ 52%]  (Sampling) 
#> Chain 4 Iteration: 1960 / 4000 [ 49%]  (Warmup) 
#> Chain 1 Iteration: 1950 / 4000 [ 48%]  (Warmup) 
#> Chain 3 Iteration: 1870 / 4000 [ 46%]  (Warmup) 
#> Chain 3 Iteration: 1880 / 4000 [ 47%]  (Warmup) 
#> Chain 2 Iteration: 2100 / 4000 [ 52%]  (Sampling) 
#> Chain 3 Iteration: 1890 / 4000 [ 47%]  (Warmup) 
#> Chain 3 Iteration: 1900 / 4000 [ 47%]  (Warmup) 
#> Chain 2 Iteration: 2110 / 4000 [ 52%]  (Sampling) 
#> Chain 3 Iteration: 1910 / 4000 [ 47%]  (Warmup) 
#> Chain 1 Iteration: 1960 / 4000 [ 49%]  (Warmup) 
#> Chain 3 Iteration: 1920 / 4000 [ 48%]  (Warmup) 
#> Chain 2 Iteration: 2120 / 4000 [ 53%]  (Sampling) 
#> Chain 4 Iteration: 1970 / 4000 [ 49%]  (Warmup) 
#> Chain 3 Iteration: 1930 / 4000 [ 48%]  (Warmup) 
#> Chain 3 Iteration: 1940 / 4000 [ 48%]  (Warmup) 
#> Chain 2 Iteration: 2130 / 4000 [ 53%]  (Sampling) 
#> Chain 3 Iteration: 1950 / 4000 [ 48%]  (Warmup) 
#> Chain 4 Iteration: 1980 / 4000 [ 49%]  (Warmup) 
#> Chain 2 Iteration: 2140 / 4000 [ 53%]  (Sampling) 
#> Chain 4 Iteration: 1990 / 4000 [ 49%]  (Warmup) 
#> Chain 1 Iteration: 1970 / 4000 [ 49%]  (Warmup) 
#> Chain 2 Iteration: 2150 / 4000 [ 53%]  (Sampling) 
#> Chain 4 Iteration: 2000 / 4000 [ 50%]  (Warmup) 
#> Chain 4 Iteration: 2001 / 4000 [ 50%]  (Sampling) 
#> Chain 3 Iteration: 1960 / 4000 [ 49%]  (Warmup) 
#> Chain 2 Iteration: 2160 / 4000 [ 54%]  (Sampling) 
#> Chain 4 Iteration: 2010 / 4000 [ 50%]  (Sampling) 
#> Chain 1 Iteration: 1980 / 4000 [ 49%]  (Warmup) 
#> Chain 2 Iteration: 2170 / 4000 [ 54%]  (Sampling) 
#> Chain 4 Iteration: 2020 / 4000 [ 50%]  (Sampling) 
#> Chain 1 Iteration: 1990 / 4000 [ 49%]  (Warmup) 
#> Chain 2 Iteration: 2180 / 4000 [ 54%]  (Sampling) 
#> Chain 4 Iteration: 2030 / 4000 [ 50%]  (Sampling) 
#> Chain 1 Iteration: 2000 / 4000 [ 50%]  (Warmup) 
#> Chain 1 Iteration: 2001 / 4000 [ 50%]  (Sampling) 
#> Chain 2 Iteration: 2190 / 4000 [ 54%]  (Sampling) 
#> Chain 4 Iteration: 2040 / 4000 [ 51%]  (Sampling) 
#> Chain 3 Iteration: 1970 / 4000 [ 49%]  (Warmup) 
#> Chain 1 Iteration: 2010 / 4000 [ 50%]  (Sampling) 
#> Chain 2 Iteration: 2200 / 4000 [ 55%]  (Sampling) 
#> Chain 4 Iteration: 2050 / 4000 [ 51%]  (Sampling) 
#> Chain 1 Iteration: 2020 / 4000 [ 50%]  (Sampling) 
#> Chain 2 Iteration: 2210 / 4000 [ 55%]  (Sampling) 
#> Chain 3 Iteration: 1980 / 4000 [ 49%]  (Warmup) 
#> Chain 4 Iteration: 2060 / 4000 [ 51%]  (Sampling) 
#> Chain 1 Iteration: 2030 / 4000 [ 50%]  (Sampling) 
#> Chain 2 Iteration: 2220 / 4000 [ 55%]  (Sampling) 
#> Chain 3 Iteration: 1990 / 4000 [ 49%]  (Warmup) 
#> Chain 4 Iteration: 2070 / 4000 [ 51%]  (Sampling) 
#> Chain 1 Iteration: 2040 / 4000 [ 51%]  (Sampling) 
#> Chain 3 Iteration: 2000 / 4000 [ 50%]  (Warmup) 
#> Chain 3 Iteration: 2001 / 4000 [ 50%]  (Sampling) 
#> Chain 2 Iteration: 2230 / 4000 [ 55%]  (Sampling) 
#> Chain 4 Iteration: 2080 / 4000 [ 52%]  (Sampling) 
#> Chain 1 Iteration: 2050 / 4000 [ 51%]  (Sampling) 
#> Chain 3 Iteration: 2010 / 4000 [ 50%]  (Sampling) 
#> Chain 2 Iteration: 2240 / 4000 [ 56%]  (Sampling) 
#> Chain 4 Iteration: 2090 / 4000 [ 52%]  (Sampling) 
#> Chain 1 Iteration: 2060 / 4000 [ 51%]  (Sampling) 
#> Chain 3 Iteration: 2020 / 4000 [ 50%]  (Sampling) 
#> Chain 2 Iteration: 2250 / 4000 [ 56%]  (Sampling) 
#> Chain 4 Iteration: 2100 / 4000 [ 52%]  (Sampling) 
#> Chain 1 Iteration: 2070 / 4000 [ 51%]  (Sampling) 
#> Chain 3 Iteration: 2030 / 4000 [ 50%]  (Sampling) 
#> Chain 2 Iteration: 2260 / 4000 [ 56%]  (Sampling) 
#> Chain 4 Iteration: 2110 / 4000 [ 52%]  (Sampling) 
#> Chain 1 Iteration: 2080 / 4000 [ 52%]  (Sampling) 
#> Chain 3 Iteration: 2040 / 4000 [ 51%]  (Sampling) 
#> Chain 2 Iteration: 2270 / 4000 [ 56%]  (Sampling) 
#> Chain 4 Iteration: 2120 / 4000 [ 53%]  (Sampling) 
#> Chain 1 Iteration: 2090 / 4000 [ 52%]  (Sampling) 
#> Chain 3 Iteration: 2050 / 4000 [ 51%]  (Sampling) 
#> Chain 2 Iteration: 2280 / 4000 [ 57%]  (Sampling) 
#> Chain 4 Iteration: 2130 / 4000 [ 53%]  (Sampling) 
#> Chain 1 Iteration: 2100 / 4000 [ 52%]  (Sampling) 
#> Chain 3 Iteration: 2060 / 4000 [ 51%]  (Sampling) 
#> Chain 2 Iteration: 2290 / 4000 [ 57%]  (Sampling) 
#> Chain 4 Iteration: 2140 / 4000 [ 53%]  (Sampling) 
#> Chain 1 Iteration: 2110 / 4000 [ 52%]  (Sampling) 
#> Chain 3 Iteration: 2070 / 4000 [ 51%]  (Sampling) 
#> Chain 2 Iteration: 2300 / 4000 [ 57%]  (Sampling) 
#> Chain 4 Iteration: 2150 / 4000 [ 53%]  (Sampling) 
#> Chain 1 Iteration: 2120 / 4000 [ 53%]  (Sampling) 
#> Chain 3 Iteration: 2080 / 4000 [ 52%]  (Sampling) 
#> Chain 2 Iteration: 2310 / 4000 [ 57%]  (Sampling) 
#> Chain 4 Iteration: 2160 / 4000 [ 54%]  (Sampling) 
#> Chain 1 Iteration: 2130 / 4000 [ 53%]  (Sampling) 
#> Chain 3 Iteration: 2090 / 4000 [ 52%]  (Sampling) 
#> Chain 2 Iteration: 2320 / 4000 [ 58%]  (Sampling) 
#> Chain 4 Iteration: 2170 / 4000 [ 54%]  (Sampling) 
#> Chain 1 Iteration: 2140 / 4000 [ 53%]  (Sampling) 
#> Chain 3 Iteration: 2100 / 4000 [ 52%]  (Sampling) 
#> Chain 2 Iteration: 2330 / 4000 [ 58%]  (Sampling) 
#> Chain 4 Iteration: 2180 / 4000 [ 54%]  (Sampling) 
#> Chain 1 Iteration: 2150 / 4000 [ 53%]  (Sampling) 
#> Chain 3 Iteration: 2110 / 4000 [ 52%]  (Sampling) 
#> Chain 2 Iteration: 2340 / 4000 [ 58%]  (Sampling) 
#> Chain 4 Iteration: 2190 / 4000 [ 54%]  (Sampling) 
#> Chain 1 Iteration: 2160 / 4000 [ 54%]  (Sampling) 
#> Chain 3 Iteration: 2120 / 4000 [ 53%]  (Sampling) 
#> Chain 2 Iteration: 2350 / 4000 [ 58%]  (Sampling) 
#> Chain 4 Iteration: 2200 / 4000 [ 55%]  (Sampling) 
#> Chain 1 Iteration: 2170 / 4000 [ 54%]  (Sampling) 
#> Chain 3 Iteration: 2130 / 4000 [ 53%]  (Sampling) 
#> Chain 2 Iteration: 2360 / 4000 [ 59%]  (Sampling) 
#> Chain 4 Iteration: 2210 / 4000 [ 55%]  (Sampling) 
#> Chain 1 Iteration: 2180 / 4000 [ 54%]  (Sampling) 
#> Chain 3 Iteration: 2140 / 4000 [ 53%]  (Sampling) 
#> Chain 2 Iteration: 2370 / 4000 [ 59%]  (Sampling) 
#> Chain 4 Iteration: 2220 / 4000 [ 55%]  (Sampling) 
#> Chain 1 Iteration: 2190 / 4000 [ 54%]  (Sampling) 
#> Chain 3 Iteration: 2150 / 4000 [ 53%]  (Sampling) 
#> Chain 2 Iteration: 2380 / 4000 [ 59%]  (Sampling) 
#> Chain 4 Iteration: 2230 / 4000 [ 55%]  (Sampling) 
#> Chain 1 Iteration: 2200 / 4000 [ 55%]  (Sampling) 
#> Chain 3 Iteration: 2160 / 4000 [ 54%]  (Sampling) 
#> Chain 2 Iteration: 2390 / 4000 [ 59%]  (Sampling) 
#> Chain 4 Iteration: 2240 / 4000 [ 56%]  (Sampling) 
#> Chain 1 Iteration: 2210 / 4000 [ 55%]  (Sampling) 
#> Chain 3 Iteration: 2170 / 4000 [ 54%]  (Sampling) 
#> Chain 2 Iteration: 2400 / 4000 [ 60%]  (Sampling) 
#> Chain 4 Iteration: 2250 / 4000 [ 56%]  (Sampling) 
#> Chain 1 Iteration: 2220 / 4000 [ 55%]  (Sampling) 
#> Chain 3 Iteration: 2180 / 4000 [ 54%]  (Sampling) 
#> Chain 2 Iteration: 2410 / 4000 [ 60%]  (Sampling) 
#> Chain 4 Iteration: 2260 / 4000 [ 56%]  (Sampling) 
#> Chain 1 Iteration: 2230 / 4000 [ 55%]  (Sampling) 
#> Chain 3 Iteration: 2190 / 4000 [ 54%]  (Sampling) 
#> Chain 2 Iteration: 2420 / 4000 [ 60%]  (Sampling) 
#> Chain 4 Iteration: 2270 / 4000 [ 56%]  (Sampling) 
#> Chain 1 Iteration: 2240 / 4000 [ 56%]  (Sampling) 
#> Chain 3 Iteration: 2200 / 4000 [ 55%]  (Sampling) 
#> Chain 2 Iteration: 2430 / 4000 [ 60%]  (Sampling) 
#> Chain 4 Iteration: 2280 / 4000 [ 57%]  (Sampling) 
#> Chain 1 Iteration: 2250 / 4000 [ 56%]  (Sampling) 
#> Chain 3 Iteration: 2210 / 4000 [ 55%]  (Sampling) 
#> Chain 2 Iteration: 2440 / 4000 [ 61%]  (Sampling) 
#> Chain 4 Iteration: 2290 / 4000 [ 57%]  (Sampling) 
#> Chain 1 Iteration: 2260 / 4000 [ 56%]  (Sampling) 
#> Chain 3 Iteration: 2220 / 4000 [ 55%]  (Sampling) 
#> Chain 2 Iteration: 2450 / 4000 [ 61%]  (Sampling) 
#> Chain 4 Iteration: 2300 / 4000 [ 57%]  (Sampling) 
#> Chain 1 Iteration: 2270 / 4000 [ 56%]  (Sampling) 
#> Chain 3 Iteration: 2230 / 4000 [ 55%]  (Sampling) 
#> Chain 2 Iteration: 2460 / 4000 [ 61%]  (Sampling) 
#> Chain 4 Iteration: 2310 / 4000 [ 57%]  (Sampling) 
#> Chain 1 Iteration: 2280 / 4000 [ 57%]  (Sampling) 
#> Chain 3 Iteration: 2240 / 4000 [ 56%]  (Sampling) 
#> Chain 2 Iteration: 2470 / 4000 [ 61%]  (Sampling) 
#> Chain 4 Iteration: 2320 / 4000 [ 58%]  (Sampling) 
#> Chain 1 Iteration: 2290 / 4000 [ 57%]  (Sampling) 
#> Chain 3 Iteration: 2250 / 4000 [ 56%]  (Sampling) 
#> Chain 2 Iteration: 2480 / 4000 [ 62%]  (Sampling) 
#> Chain 4 Iteration: 2330 / 4000 [ 58%]  (Sampling) 
#> Chain 1 Iteration: 2300 / 4000 [ 57%]  (Sampling) 
#> Chain 3 Iteration: 2260 / 4000 [ 56%]  (Sampling) 
#> Chain 2 Iteration: 2490 / 4000 [ 62%]  (Sampling) 
#> Chain 4 Iteration: 2340 / 4000 [ 58%]  (Sampling) 
#> Chain 1 Iteration: 2310 / 4000 [ 57%]  (Sampling) 
#> Chain 3 Iteration: 2270 / 4000 [ 56%]  (Sampling) 
#> Chain 2 Iteration: 2500 / 4000 [ 62%]  (Sampling) 
#> Chain 4 Iteration: 2350 / 4000 [ 58%]  (Sampling) 
#> Chain 1 Iteration: 2320 / 4000 [ 58%]  (Sampling) 
#> Chain 3 Iteration: 2280 / 4000 [ 57%]  (Sampling) 
#> Chain 2 Iteration: 2510 / 4000 [ 62%]  (Sampling) 
#> Chain 4 Iteration: 2360 / 4000 [ 59%]  (Sampling) 
#> Chain 1 Iteration: 2330 / 4000 [ 58%]  (Sampling) 
#> Chain 3 Iteration: 2290 / 4000 [ 57%]  (Sampling) 
#> Chain 2 Iteration: 2520 / 4000 [ 63%]  (Sampling) 
#> Chain 4 Iteration: 2370 / 4000 [ 59%]  (Sampling) 
#> Chain 1 Iteration: 2340 / 4000 [ 58%]  (Sampling) 
#> Chain 3 Iteration: 2300 / 4000 [ 57%]  (Sampling) 
#> Chain 2 Iteration: 2530 / 4000 [ 63%]  (Sampling) 
#> Chain 4 Iteration: 2380 / 4000 [ 59%]  (Sampling) 
#> Chain 1 Iteration: 2350 / 4000 [ 58%]  (Sampling) 
#> Chain 3 Iteration: 2310 / 4000 [ 57%]  (Sampling) 
#> Chain 2 Iteration: 2540 / 4000 [ 63%]  (Sampling) 
#> Chain 4 Iteration: 2390 / 4000 [ 59%]  (Sampling) 
#> Chain 1 Iteration: 2360 / 4000 [ 59%]  (Sampling) 
#> Chain 3 Iteration: 2320 / 4000 [ 58%]  (Sampling) 
#> Chain 2 Iteration: 2550 / 4000 [ 63%]  (Sampling) 
#> Chain 4 Iteration: 2400 / 4000 [ 60%]  (Sampling) 
#> Chain 1 Iteration: 2370 / 4000 [ 59%]  (Sampling) 
#> Chain 3 Iteration: 2330 / 4000 [ 58%]  (Sampling) 
#> Chain 2 Iteration: 2560 / 4000 [ 64%]  (Sampling) 
#> Chain 4 Iteration: 2410 / 4000 [ 60%]  (Sampling) 
#> Chain 1 Iteration: 2380 / 4000 [ 59%]  (Sampling) 
#> Chain 3 Iteration: 2340 / 4000 [ 58%]  (Sampling) 
#> Chain 2 Iteration: 2570 / 4000 [ 64%]  (Sampling) 
#> Chain 4 Iteration: 2420 / 4000 [ 60%]  (Sampling) 
#> Chain 1 Iteration: 2390 / 4000 [ 59%]  (Sampling) 
#> Chain 3 Iteration: 2350 / 4000 [ 58%]  (Sampling) 
#> Chain 2 Iteration: 2580 / 4000 [ 64%]  (Sampling) 
#> Chain 4 Iteration: 2430 / 4000 [ 60%]  (Sampling) 
#> Chain 1 Iteration: 2400 / 4000 [ 60%]  (Sampling) 
#> Chain 3 Iteration: 2360 / 4000 [ 59%]  (Sampling) 
#> Chain 2 Iteration: 2590 / 4000 [ 64%]  (Sampling) 
#> Chain 4 Iteration: 2440 / 4000 [ 61%]  (Sampling) 
#> Chain 1 Iteration: 2410 / 4000 [ 60%]  (Sampling) 
#> Chain 3 Iteration: 2370 / 4000 [ 59%]  (Sampling) 
#> Chain 2 Iteration: 2600 / 4000 [ 65%]  (Sampling) 
#> Chain 4 Iteration: 2450 / 4000 [ 61%]  (Sampling) 
#> Chain 1 Iteration: 2420 / 4000 [ 60%]  (Sampling) 
#> Chain 3 Iteration: 2380 / 4000 [ 59%]  (Sampling) 
#> Chain 2 Iteration: 2610 / 4000 [ 65%]  (Sampling) 
#> Chain 4 Iteration: 2460 / 4000 [ 61%]  (Sampling) 
#> Chain 1 Iteration: 2430 / 4000 [ 60%]  (Sampling) 
#> Chain 3 Iteration: 2390 / 4000 [ 59%]  (Sampling) 
#> Chain 2 Iteration: 2620 / 4000 [ 65%]  (Sampling) 
#> Chain 4 Iteration: 2470 / 4000 [ 61%]  (Sampling) 
#> Chain 1 Iteration: 2440 / 4000 [ 61%]  (Sampling) 
#> Chain 3 Iteration: 2400 / 4000 [ 60%]  (Sampling) 
#> Chain 2 Iteration: 2630 / 4000 [ 65%]  (Sampling) 
#> Chain 4 Iteration: 2480 / 4000 [ 62%]  (Sampling) 
#> Chain 1 Iteration: 2450 / 4000 [ 61%]  (Sampling) 
#> Chain 3 Iteration: 2410 / 4000 [ 60%]  (Sampling) 
#> Chain 2 Iteration: 2640 / 4000 [ 66%]  (Sampling) 
#> Chain 4 Iteration: 2490 / 4000 [ 62%]  (Sampling) 
#> Chain 1 Iteration: 2460 / 4000 [ 61%]  (Sampling) 
#> Chain 3 Iteration: 2420 / 4000 [ 60%]  (Sampling) 
#> Chain 2 Iteration: 2650 / 4000 [ 66%]  (Sampling) 
#> Chain 4 Iteration: 2500 / 4000 [ 62%]  (Sampling) 
#> Chain 1 Iteration: 2470 / 4000 [ 61%]  (Sampling) 
#> Chain 3 Iteration: 2430 / 4000 [ 60%]  (Sampling) 
#> Chain 2 Iteration: 2660 / 4000 [ 66%]  (Sampling) 
#> Chain 4 Iteration: 2510 / 4000 [ 62%]  (Sampling) 
#> Chain 1 Iteration: 2480 / 4000 [ 62%]  (Sampling) 
#> Chain 3 Iteration: 2440 / 4000 [ 61%]  (Sampling) 
#> Chain 2 Iteration: 2670 / 4000 [ 66%]  (Sampling) 
#> Chain 4 Iteration: 2520 / 4000 [ 63%]  (Sampling) 
#> Chain 1 Iteration: 2490 / 4000 [ 62%]  (Sampling) 
#> Chain 3 Iteration: 2450 / 4000 [ 61%]  (Sampling) 
#> Chain 2 Iteration: 2680 / 4000 [ 67%]  (Sampling) 
#> Chain 4 Iteration: 2530 / 4000 [ 63%]  (Sampling) 
#> Chain 1 Iteration: 2500 / 4000 [ 62%]  (Sampling) 
#> Chain 3 Iteration: 2460 / 4000 [ 61%]  (Sampling) 
#> Chain 2 Iteration: 2690 / 4000 [ 67%]  (Sampling) 
#> Chain 4 Iteration: 2540 / 4000 [ 63%]  (Sampling) 
#> Chain 1 Iteration: 2510 / 4000 [ 62%]  (Sampling) 
#> Chain 3 Iteration: 2470 / 4000 [ 61%]  (Sampling) 
#> Chain 2 Iteration: 2700 / 4000 [ 67%]  (Sampling) 
#> Chain 4 Iteration: 2550 / 4000 [ 63%]  (Sampling) 
#> Chain 1 Iteration: 2520 / 4000 [ 63%]  (Sampling) 
#> Chain 3 Iteration: 2480 / 4000 [ 62%]  (Sampling) 
#> Chain 2 Iteration: 2710 / 4000 [ 67%]  (Sampling) 
#> Chain 4 Iteration: 2560 / 4000 [ 64%]  (Sampling) 
#> Chain 1 Iteration: 2530 / 4000 [ 63%]  (Sampling) 
#> Chain 3 Iteration: 2490 / 4000 [ 62%]  (Sampling) 
#> Chain 2 Iteration: 2720 / 4000 [ 68%]  (Sampling) 
#> Chain 4 Iteration: 2570 / 4000 [ 64%]  (Sampling) 
#> Chain 1 Iteration: 2540 / 4000 [ 63%]  (Sampling) 
#> Chain 3 Iteration: 2500 / 4000 [ 62%]  (Sampling) 
#> Chain 2 Iteration: 2730 / 4000 [ 68%]  (Sampling) 
#> Chain 4 Iteration: 2580 / 4000 [ 64%]  (Sampling) 
#> Chain 1 Iteration: 2550 / 4000 [ 63%]  (Sampling) 
#> Chain 3 Iteration: 2510 / 4000 [ 62%]  (Sampling) 
#> Chain 2 Iteration: 2740 / 4000 [ 68%]  (Sampling) 
#> Chain 4 Iteration: 2590 / 4000 [ 64%]  (Sampling) 
#> Chain 1 Iteration: 2560 / 4000 [ 64%]  (Sampling) 
#> Chain 3 Iteration: 2520 / 4000 [ 63%]  (Sampling) 
#> Chain 2 Iteration: 2750 / 4000 [ 68%]  (Sampling) 
#> Chain 4 Iteration: 2600 / 4000 [ 65%]  (Sampling) 
#> Chain 1 Iteration: 2570 / 4000 [ 64%]  (Sampling) 
#> Chain 3 Iteration: 2530 / 4000 [ 63%]  (Sampling) 
#> Chain 2 Iteration: 2760 / 4000 [ 69%]  (Sampling) 
#> Chain 4 Iteration: 2610 / 4000 [ 65%]  (Sampling) 
#> Chain 1 Iteration: 2580 / 4000 [ 64%]  (Sampling) 
#> Chain 3 Iteration: 2540 / 4000 [ 63%]  (Sampling) 
#> Chain 2 Iteration: 2770 / 4000 [ 69%]  (Sampling) 
#> Chain 4 Iteration: 2620 / 4000 [ 65%]  (Sampling) 
#> Chain 1 Iteration: 2590 / 4000 [ 64%]  (Sampling) 
#> Chain 3 Iteration: 2550 / 4000 [ 63%]  (Sampling) 
#> Chain 2 Iteration: 2780 / 4000 [ 69%]  (Sampling) 
#> Chain 4 Iteration: 2630 / 4000 [ 65%]  (Sampling) 
#> Chain 1 Iteration: 2600 / 4000 [ 65%]  (Sampling) 
#> Chain 3 Iteration: 2560 / 4000 [ 64%]  (Sampling) 
#> Chain 2 Iteration: 2790 / 4000 [ 69%]  (Sampling) 
#> Chain 4 Iteration: 2640 / 4000 [ 66%]  (Sampling) 
#> Chain 1 Iteration: 2610 / 4000 [ 65%]  (Sampling) 
#> Chain 3 Iteration: 2570 / 4000 [ 64%]  (Sampling) 
#> Chain 2 Iteration: 2800 / 4000 [ 70%]  (Sampling) 
#> Chain 4 Iteration: 2650 / 4000 [ 66%]  (Sampling) 
#> Chain 1 Iteration: 2620 / 4000 [ 65%]  (Sampling) 
#> Chain 3 Iteration: 2580 / 4000 [ 64%]  (Sampling) 
#> Chain 2 Iteration: 2810 / 4000 [ 70%]  (Sampling) 
#> Chain 4 Iteration: 2660 / 4000 [ 66%]  (Sampling) 
#> Chain 1 Iteration: 2630 / 4000 [ 65%]  (Sampling) 
#> Chain 3 Iteration: 2590 / 4000 [ 64%]  (Sampling) 
#> Chain 2 Iteration: 2820 / 4000 [ 70%]  (Sampling) 
#> Chain 4 Iteration: 2670 / 4000 [ 66%]  (Sampling) 
#> Chain 1 Iteration: 2640 / 4000 [ 66%]  (Sampling) 
#> Chain 3 Iteration: 2600 / 4000 [ 65%]  (Sampling) 
#> Chain 2 Iteration: 2830 / 4000 [ 70%]  (Sampling) 
#> Chain 4 Iteration: 2680 / 4000 [ 67%]  (Sampling) 
#> Chain 1 Iteration: 2650 / 4000 [ 66%]  (Sampling) 
#> Chain 3 Iteration: 2610 / 4000 [ 65%]  (Sampling) 
#> Chain 2 Iteration: 2840 / 4000 [ 71%]  (Sampling) 
#> Chain 4 Iteration: 2690 / 4000 [ 67%]  (Sampling) 
#> Chain 1 Iteration: 2660 / 4000 [ 66%]  (Sampling) 
#> Chain 3 Iteration: 2620 / 4000 [ 65%]  (Sampling) 
#> Chain 2 Iteration: 2850 / 4000 [ 71%]  (Sampling) 
#> Chain 4 Iteration: 2700 / 4000 [ 67%]  (Sampling) 
#> Chain 1 Iteration: 2670 / 4000 [ 66%]  (Sampling) 
#> Chain 3 Iteration: 2630 / 4000 [ 65%]  (Sampling) 
#> Chain 2 Iteration: 2860 / 4000 [ 71%]  (Sampling) 
#> Chain 4 Iteration: 2710 / 4000 [ 67%]  (Sampling) 
#> Chain 1 Iteration: 2680 / 4000 [ 67%]  (Sampling) 
#> Chain 3 Iteration: 2640 / 4000 [ 66%]  (Sampling) 
#> Chain 2 Iteration: 2870 / 4000 [ 71%]  (Sampling) 
#> Chain 4 Iteration: 2720 / 4000 [ 68%]  (Sampling) 
#> Chain 1 Iteration: 2690 / 4000 [ 67%]  (Sampling) 
#> Chain 3 Iteration: 2650 / 4000 [ 66%]  (Sampling) 
#> Chain 2 Iteration: 2880 / 4000 [ 72%]  (Sampling) 
#> Chain 4 Iteration: 2730 / 4000 [ 68%]  (Sampling) 
#> Chain 1 Iteration: 2700 / 4000 [ 67%]  (Sampling) 
#> Chain 3 Iteration: 2660 / 4000 [ 66%]  (Sampling) 
#> Chain 2 Iteration: 2890 / 4000 [ 72%]  (Sampling) 
#> Chain 4 Iteration: 2740 / 4000 [ 68%]  (Sampling) 
#> Chain 1 Iteration: 2710 / 4000 [ 67%]  (Sampling) 
#> Chain 3 Iteration: 2670 / 4000 [ 66%]  (Sampling) 
#> Chain 2 Iteration: 2900 / 4000 [ 72%]  (Sampling) 
#> Chain 4 Iteration: 2750 / 4000 [ 68%]  (Sampling) 
#> Chain 1 Iteration: 2720 / 4000 [ 68%]  (Sampling) 
#> Chain 3 Iteration: 2680 / 4000 [ 67%]  (Sampling) 
#> Chain 2 Iteration: 2910 / 4000 [ 72%]  (Sampling) 
#> Chain 4 Iteration: 2760 / 4000 [ 69%]  (Sampling) 
#> Chain 1 Iteration: 2730 / 4000 [ 68%]  (Sampling) 
#> Chain 3 Iteration: 2690 / 4000 [ 67%]  (Sampling) 
#> Chain 2 Iteration: 2920 / 4000 [ 73%]  (Sampling) 
#> Chain 4 Iteration: 2770 / 4000 [ 69%]  (Sampling) 
#> Chain 1 Iteration: 2740 / 4000 [ 68%]  (Sampling) 
#> Chain 3 Iteration: 2700 / 4000 [ 67%]  (Sampling) 
#> Chain 2 Iteration: 2930 / 4000 [ 73%]  (Sampling) 
#> Chain 4 Iteration: 2780 / 4000 [ 69%]  (Sampling) 
#> Chain 1 Iteration: 2750 / 4000 [ 68%]  (Sampling) 
#> Chain 3 Iteration: 2710 / 4000 [ 67%]  (Sampling) 
#> Chain 2 Iteration: 2940 / 4000 [ 73%]  (Sampling) 
#> Chain 4 Iteration: 2790 / 4000 [ 69%]  (Sampling) 
#> Chain 1 Iteration: 2760 / 4000 [ 69%]  (Sampling) 
#> Chain 3 Iteration: 2720 / 4000 [ 68%]  (Sampling) 
#> Chain 2 Iteration: 2950 / 4000 [ 73%]  (Sampling) 
#> Chain 4 Iteration: 2800 / 4000 [ 70%]  (Sampling) 
#> Chain 1 Iteration: 2770 / 4000 [ 69%]  (Sampling) 
#> Chain 3 Iteration: 2730 / 4000 [ 68%]  (Sampling) 
#> Chain 2 Iteration: 2960 / 4000 [ 74%]  (Sampling) 
#> Chain 4 Iteration: 2810 / 4000 [ 70%]  (Sampling) 
#> Chain 1 Iteration: 2780 / 4000 [ 69%]  (Sampling) 
#> Chain 3 Iteration: 2740 / 4000 [ 68%]  (Sampling) 
#> Chain 2 Iteration: 2970 / 4000 [ 74%]  (Sampling) 
#> Chain 4 Iteration: 2820 / 4000 [ 70%]  (Sampling) 
#> Chain 1 Iteration: 2790 / 4000 [ 69%]  (Sampling) 
#> Chain 3 Iteration: 2750 / 4000 [ 68%]  (Sampling) 
#> Chain 2 Iteration: 2980 / 4000 [ 74%]  (Sampling) 
#> Chain 4 Iteration: 2830 / 4000 [ 70%]  (Sampling) 
#> Chain 1 Iteration: 2800 / 4000 [ 70%]  (Sampling) 
#> Chain 3 Iteration: 2760 / 4000 [ 69%]  (Sampling) 
#> Chain 2 Iteration: 2990 / 4000 [ 74%]  (Sampling) 
#> Chain 4 Iteration: 2840 / 4000 [ 71%]  (Sampling) 
#> Chain 1 Iteration: 2810 / 4000 [ 70%]  (Sampling) 
#> Chain 3 Iteration: 2770 / 4000 [ 69%]  (Sampling) 
#> Chain 2 Iteration: 3000 / 4000 [ 75%]  (Sampling) 
#> Chain 4 Iteration: 2850 / 4000 [ 71%]  (Sampling) 
#> Chain 1 Iteration: 2820 / 4000 [ 70%]  (Sampling) 
#> Chain 3 Iteration: 2780 / 4000 [ 69%]  (Sampling) 
#> Chain 2 Iteration: 3010 / 4000 [ 75%]  (Sampling) 
#> Chain 4 Iteration: 2860 / 4000 [ 71%]  (Sampling) 
#> Chain 1 Iteration: 2830 / 4000 [ 70%]  (Sampling) 
#> Chain 3 Iteration: 2790 / 4000 [ 69%]  (Sampling) 
#> Chain 2 Iteration: 3020 / 4000 [ 75%]  (Sampling) 
#> Chain 4 Iteration: 2870 / 4000 [ 71%]  (Sampling) 
#> Chain 1 Iteration: 2840 / 4000 [ 71%]  (Sampling) 
#> Chain 3 Iteration: 2800 / 4000 [ 70%]  (Sampling) 
#> Chain 2 Iteration: 3030 / 4000 [ 75%]  (Sampling) 
#> Chain 4 Iteration: 2880 / 4000 [ 72%]  (Sampling) 
#> Chain 1 Iteration: 2850 / 4000 [ 71%]  (Sampling) 
#> Chain 3 Iteration: 2810 / 4000 [ 70%]  (Sampling) 
#> Chain 2 Iteration: 3040 / 4000 [ 76%]  (Sampling) 
#> Chain 4 Iteration: 2890 / 4000 [ 72%]  (Sampling) 
#> Chain 1 Iteration: 2860 / 4000 [ 71%]  (Sampling) 
#> Chain 3 Iteration: 2820 / 4000 [ 70%]  (Sampling) 
#> Chain 2 Iteration: 3050 / 4000 [ 76%]  (Sampling) 
#> Chain 4 Iteration: 2900 / 4000 [ 72%]  (Sampling) 
#> Chain 1 Iteration: 2870 / 4000 [ 71%]  (Sampling) 
#> Chain 3 Iteration: 2830 / 4000 [ 70%]  (Sampling) 
#> Chain 2 Iteration: 3060 / 4000 [ 76%]  (Sampling) 
#> Chain 4 Iteration: 2910 / 4000 [ 72%]  (Sampling) 
#> Chain 1 Iteration: 2880 / 4000 [ 72%]  (Sampling) 
#> Chain 3 Iteration: 2840 / 4000 [ 71%]  (Sampling) 
#> Chain 2 Iteration: 3070 / 4000 [ 76%]  (Sampling) 
#> Chain 4 Iteration: 2920 / 4000 [ 73%]  (Sampling) 
#> Chain 1 Iteration: 2890 / 4000 [ 72%]  (Sampling) 
#> Chain 3 Iteration: 2850 / 4000 [ 71%]  (Sampling) 
#> Chain 2 Iteration: 3080 / 4000 [ 77%]  (Sampling) 
#> Chain 4 Iteration: 2930 / 4000 [ 73%]  (Sampling) 
#> Chain 1 Iteration: 2900 / 4000 [ 72%]  (Sampling) 
#> Chain 3 Iteration: 2860 / 4000 [ 71%]  (Sampling) 
#> Chain 2 Iteration: 3090 / 4000 [ 77%]  (Sampling) 
#> Chain 4 Iteration: 2940 / 4000 [ 73%]  (Sampling) 
#> Chain 1 Iteration: 2910 / 4000 [ 72%]  (Sampling) 
#> Chain 3 Iteration: 2870 / 4000 [ 71%]  (Sampling) 
#> Chain 2 Iteration: 3100 / 4000 [ 77%]  (Sampling) 
#> Chain 4 Iteration: 2950 / 4000 [ 73%]  (Sampling) 
#> Chain 1 Iteration: 2920 / 4000 [ 73%]  (Sampling) 
#> Chain 3 Iteration: 2880 / 4000 [ 72%]  (Sampling) 
#> Chain 2 Iteration: 3110 / 4000 [ 77%]  (Sampling) 
#> Chain 4 Iteration: 2960 / 4000 [ 74%]  (Sampling) 
#> Chain 1 Iteration: 2930 / 4000 [ 73%]  (Sampling) 
#> Chain 3 Iteration: 2890 / 4000 [ 72%]  (Sampling) 
#> Chain 2 Iteration: 3120 / 4000 [ 78%]  (Sampling) 
#> Chain 4 Iteration: 2970 / 4000 [ 74%]  (Sampling) 
#> Chain 1 Iteration: 2940 / 4000 [ 73%]  (Sampling) 
#> Chain 3 Iteration: 2900 / 4000 [ 72%]  (Sampling) 
#> Chain 2 Iteration: 3130 / 4000 [ 78%]  (Sampling) 
#> Chain 4 Iteration: 2980 / 4000 [ 74%]  (Sampling) 
#> Chain 1 Iteration: 2950 / 4000 [ 73%]  (Sampling) 
#> Chain 3 Iteration: 2910 / 4000 [ 72%]  (Sampling) 
#> Chain 2 Iteration: 3140 / 4000 [ 78%]  (Sampling) 
#> Chain 4 Iteration: 2990 / 4000 [ 74%]  (Sampling) 
#> Chain 1 Iteration: 2960 / 4000 [ 74%]  (Sampling) 
#> Chain 3 Iteration: 2920 / 4000 [ 73%]  (Sampling) 
#> Chain 2 Iteration: 3150 / 4000 [ 78%]  (Sampling) 
#> Chain 4 Iteration: 3000 / 4000 [ 75%]  (Sampling) 
#> Chain 1 Iteration: 2970 / 4000 [ 74%]  (Sampling) 
#> Chain 3 Iteration: 2930 / 4000 [ 73%]  (Sampling) 
#> Chain 2 Iteration: 3160 / 4000 [ 79%]  (Sampling) 
#> Chain 4 Iteration: 3010 / 4000 [ 75%]  (Sampling) 
#> Chain 1 Iteration: 2980 / 4000 [ 74%]  (Sampling) 
#> Chain 3 Iteration: 2940 / 4000 [ 73%]  (Sampling) 
#> Chain 2 Iteration: 3170 / 4000 [ 79%]  (Sampling) 
#> Chain 4 Iteration: 3020 / 4000 [ 75%]  (Sampling) 
#> Chain 1 Iteration: 2990 / 4000 [ 74%]  (Sampling) 
#> Chain 3 Iteration: 2950 / 4000 [ 73%]  (Sampling) 
#> Chain 2 Iteration: 3180 / 4000 [ 79%]  (Sampling) 
#> Chain 4 Iteration: 3030 / 4000 [ 75%]  (Sampling) 
#> Chain 1 Iteration: 3000 / 4000 [ 75%]  (Sampling) 
#> Chain 3 Iteration: 2960 / 4000 [ 74%]  (Sampling) 
#> Chain 2 Iteration: 3190 / 4000 [ 79%]  (Sampling) 
#> Chain 4 Iteration: 3040 / 4000 [ 76%]  (Sampling) 
#> Chain 1 Iteration: 3010 / 4000 [ 75%]  (Sampling) 
#> Chain 3 Iteration: 2970 / 4000 [ 74%]  (Sampling) 
#> Chain 2 Iteration: 3200 / 4000 [ 80%]  (Sampling) 
#> Chain 4 Iteration: 3050 / 4000 [ 76%]  (Sampling) 
#> Chain 1 Iteration: 3020 / 4000 [ 75%]  (Sampling) 
#> Chain 2 Iteration: 3210 / 4000 [ 80%]  (Sampling) 
#> Chain 3 Iteration: 2980 / 4000 [ 74%]  (Sampling) 
#> Chain 4 Iteration: 3060 / 4000 [ 76%]  (Sampling) 
#> Chain 1 Iteration: 3030 / 4000 [ 75%]  (Sampling) 
#> Chain 2 Iteration: 3220 / 4000 [ 80%]  (Sampling) 
#> Chain 3 Iteration: 2990 / 4000 [ 74%]  (Sampling) 
#> Chain 4 Iteration: 3070 / 4000 [ 76%]  (Sampling) 
#> Chain 1 Iteration: 3040 / 4000 [ 76%]  (Sampling) 
#> Chain 2 Iteration: 3230 / 4000 [ 80%]  (Sampling) 
#> Chain 3 Iteration: 3000 / 4000 [ 75%]  (Sampling) 
#> Chain 4 Iteration: 3080 / 4000 [ 77%]  (Sampling) 
#> Chain 1 Iteration: 3050 / 4000 [ 76%]  (Sampling) 
#> Chain 2 Iteration: 3240 / 4000 [ 81%]  (Sampling) 
#> Chain 3 Iteration: 3010 / 4000 [ 75%]  (Sampling) 
#> Chain 4 Iteration: 3090 / 4000 [ 77%]  (Sampling) 
#> Chain 1 Iteration: 3060 / 4000 [ 76%]  (Sampling) 
#> Chain 2 Iteration: 3250 / 4000 [ 81%]  (Sampling) 
#> Chain 3 Iteration: 3020 / 4000 [ 75%]  (Sampling) 
#> Chain 4 Iteration: 3100 / 4000 [ 77%]  (Sampling) 
#> Chain 1 Iteration: 3070 / 4000 [ 76%]  (Sampling) 
#> Chain 2 Iteration: 3260 / 4000 [ 81%]  (Sampling) 
#> Chain 3 Iteration: 3030 / 4000 [ 75%]  (Sampling) 
#> Chain 4 Iteration: 3110 / 4000 [ 77%]  (Sampling) 
#> Chain 1 Iteration: 3080 / 4000 [ 77%]  (Sampling) 
#> Chain 2 Iteration: 3270 / 4000 [ 81%]  (Sampling) 
#> Chain 3 Iteration: 3040 / 4000 [ 76%]  (Sampling) 
#> Chain 4 Iteration: 3120 / 4000 [ 78%]  (Sampling) 
#> Chain 1 Iteration: 3090 / 4000 [ 77%]  (Sampling) 
#> Chain 2 Iteration: 3280 / 4000 [ 82%]  (Sampling) 
#> Chain 3 Iteration: 3050 / 4000 [ 76%]  (Sampling) 
#> Chain 4 Iteration: 3130 / 4000 [ 78%]  (Sampling) 
#> Chain 1 Iteration: 3100 / 4000 [ 77%]  (Sampling) 
#> Chain 2 Iteration: 3290 / 4000 [ 82%]  (Sampling) 
#> Chain 3 Iteration: 3060 / 4000 [ 76%]  (Sampling) 
#> Chain 4 Iteration: 3140 / 4000 [ 78%]  (Sampling) 
#> Chain 1 Iteration: 3110 / 4000 [ 77%]  (Sampling) 
#> Chain 2 Iteration: 3300 / 4000 [ 82%]  (Sampling) 
#> Chain 3 Iteration: 3070 / 4000 [ 76%]  (Sampling) 
#> Chain 4 Iteration: 3150 / 4000 [ 78%]  (Sampling) 
#> Chain 1 Iteration: 3120 / 4000 [ 78%]  (Sampling) 
#> Chain 2 Iteration: 3310 / 4000 [ 82%]  (Sampling) 
#> Chain 3 Iteration: 3080 / 4000 [ 77%]  (Sampling) 
#> Chain 4 Iteration: 3160 / 4000 [ 79%]  (Sampling) 
#> Chain 1 Iteration: 3130 / 4000 [ 78%]  (Sampling) 
#> Chain 2 Iteration: 3320 / 4000 [ 83%]  (Sampling) 
#> Chain 3 Iteration: 3090 / 4000 [ 77%]  (Sampling) 
#> Chain 4 Iteration: 3170 / 4000 [ 79%]  (Sampling) 
#> Chain 1 Iteration: 3140 / 4000 [ 78%]  (Sampling) 
#> Chain 2 Iteration: 3330 / 4000 [ 83%]  (Sampling) 
#> Chain 3 Iteration: 3100 / 4000 [ 77%]  (Sampling) 
#> Chain 4 Iteration: 3180 / 4000 [ 79%]  (Sampling) 
#> Chain 1 Iteration: 3150 / 4000 [ 78%]  (Sampling) 
#> Chain 2 Iteration: 3340 / 4000 [ 83%]  (Sampling) 
#> Chain 3 Iteration: 3110 / 4000 [ 77%]  (Sampling) 
#> Chain 4 Iteration: 3190 / 4000 [ 79%]  (Sampling) 
#> Chain 1 Iteration: 3160 / 4000 [ 79%]  (Sampling) 
#> Chain 2 Iteration: 3350 / 4000 [ 83%]  (Sampling) 
#> Chain 3 Iteration: 3120 / 4000 [ 78%]  (Sampling) 
#> Chain 4 Iteration: 3200 / 4000 [ 80%]  (Sampling) 
#> Chain 1 Iteration: 3170 / 4000 [ 79%]  (Sampling) 
#> Chain 2 Iteration: 3360 / 4000 [ 84%]  (Sampling) 
#> Chain 3 Iteration: 3130 / 4000 [ 78%]  (Sampling) 
#> Chain 4 Iteration: 3210 / 4000 [ 80%]  (Sampling) 
#> Chain 1 Iteration: 3180 / 4000 [ 79%]  (Sampling) 
#> Chain 2 Iteration: 3370 / 4000 [ 84%]  (Sampling) 
#> Chain 3 Iteration: 3140 / 4000 [ 78%]  (Sampling) 
#> Chain 4 Iteration: 3220 / 4000 [ 80%]  (Sampling) 
#> Chain 1 Iteration: 3190 / 4000 [ 79%]  (Sampling) 
#> Chain 2 Iteration: 3380 / 4000 [ 84%]  (Sampling) 
#> Chain 3 Iteration: 3150 / 4000 [ 78%]  (Sampling) 
#> Chain 4 Iteration: 3230 / 4000 [ 80%]  (Sampling) 
#> Chain 1 Iteration: 3200 / 4000 [ 80%]  (Sampling) 
#> Chain 2 Iteration: 3390 / 4000 [ 84%]  (Sampling) 
#> Chain 3 Iteration: 3160 / 4000 [ 79%]  (Sampling) 
#> Chain 4 Iteration: 3240 / 4000 [ 81%]  (Sampling) 
#> Chain 1 Iteration: 3210 / 4000 [ 80%]  (Sampling) 
#> Chain 2 Iteration: 3400 / 4000 [ 85%]  (Sampling) 
#> Chain 3 Iteration: 3170 / 4000 [ 79%]  (Sampling) 
#> Chain 4 Iteration: 3250 / 4000 [ 81%]  (Sampling) 
#> Chain 1 Iteration: 3220 / 4000 [ 80%]  (Sampling) 
#> Chain 3 Iteration: 3180 / 4000 [ 79%]  (Sampling) 
#> Chain 2 Iteration: 3410 / 4000 [ 85%]  (Sampling) 
#> Chain 4 Iteration: 3260 / 4000 [ 81%]  (Sampling) 
#> Chain 1 Iteration: 3230 / 4000 [ 80%]  (Sampling) 
#> Chain 2 Iteration: 3420 / 4000 [ 85%]  (Sampling) 
#> Chain 3 Iteration: 3190 / 4000 [ 79%]  (Sampling) 
#> Chain 4 Iteration: 3270 / 4000 [ 81%]  (Sampling) 
#> Chain 1 Iteration: 3240 / 4000 [ 81%]  (Sampling) 
#> Chain 2 Iteration: 3430 / 4000 [ 85%]  (Sampling) 
#> Chain 3 Iteration: 3200 / 4000 [ 80%]  (Sampling) 
#> Chain 4 Iteration: 3280 / 4000 [ 82%]  (Sampling) 
#> Chain 1 Iteration: 3250 / 4000 [ 81%]  (Sampling) 
#> Chain 3 Iteration: 3210 / 4000 [ 80%]  (Sampling) 
#> Chain 2 Iteration: 3440 / 4000 [ 86%]  (Sampling) 
#> Chain 4 Iteration: 3290 / 4000 [ 82%]  (Sampling) 
#> Chain 1 Iteration: 3260 / 4000 [ 81%]  (Sampling) 
#> Chain 3 Iteration: 3220 / 4000 [ 80%]  (Sampling) 
#> Chain 2 Iteration: 3450 / 4000 [ 86%]  (Sampling) 
#> Chain 4 Iteration: 3300 / 4000 [ 82%]  (Sampling) 
#> Chain 1 Iteration: 3270 / 4000 [ 81%]  (Sampling) 
#> Chain 3 Iteration: 3230 / 4000 [ 80%]  (Sampling) 
#> Chain 2 Iteration: 3460 / 4000 [ 86%]  (Sampling) 
#> Chain 4 Iteration: 3310 / 4000 [ 82%]  (Sampling) 
#> Chain 1 Iteration: 3280 / 4000 [ 82%]  (Sampling) 
#> Chain 3 Iteration: 3240 / 4000 [ 81%]  (Sampling) 
#> Chain 2 Iteration: 3470 / 4000 [ 86%]  (Sampling) 
#> Chain 4 Iteration: 3320 / 4000 [ 83%]  (Sampling) 
#> Chain 1 Iteration: 3290 / 4000 [ 82%]  (Sampling) 
#> Chain 3 Iteration: 3250 / 4000 [ 81%]  (Sampling) 
#> Chain 2 Iteration: 3480 / 4000 [ 87%]  (Sampling) 
#> Chain 4 Iteration: 3330 / 4000 [ 83%]  (Sampling) 
#> Chain 1 Iteration: 3300 / 4000 [ 82%]  (Sampling) 
#> Chain 3 Iteration: 3260 / 4000 [ 81%]  (Sampling) 
#> Chain 2 Iteration: 3490 / 4000 [ 87%]  (Sampling) 
#> Chain 4 Iteration: 3340 / 4000 [ 83%]  (Sampling) 
#> Chain 1 Iteration: 3310 / 4000 [ 82%]  (Sampling) 
#> Chain 3 Iteration: 3270 / 4000 [ 81%]  (Sampling) 
#> Chain 2 Iteration: 3500 / 4000 [ 87%]  (Sampling) 
#> Chain 4 Iteration: 3350 / 4000 [ 83%]  (Sampling) 
#> Chain 1 Iteration: 3320 / 4000 [ 83%]  (Sampling) 
#> Chain 3 Iteration: 3280 / 4000 [ 82%]  (Sampling) 
#> Chain 2 Iteration: 3510 / 4000 [ 87%]  (Sampling) 
#> Chain 4 Iteration: 3360 / 4000 [ 84%]  (Sampling) 
#> Chain 1 Iteration: 3330 / 4000 [ 83%]  (Sampling) 
#> Chain 2 Iteration: 3520 / 4000 [ 88%]  (Sampling) 
#> Chain 3 Iteration: 3290 / 4000 [ 82%]  (Sampling) 
#> Chain 4 Iteration: 3370 / 4000 [ 84%]  (Sampling) 
#> Chain 1 Iteration: 3340 / 4000 [ 83%]  (Sampling) 
#> Chain 3 Iteration: 3300 / 4000 [ 82%]  (Sampling) 
#> Chain 2 Iteration: 3530 / 4000 [ 88%]  (Sampling) 
#> Chain 4 Iteration: 3380 / 4000 [ 84%]  (Sampling) 
#> Chain 1 Iteration: 3350 / 4000 [ 83%]  (Sampling) 
#> Chain 2 Iteration: 3540 / 4000 [ 88%]  (Sampling) 
#> Chain 3 Iteration: 3310 / 4000 [ 82%]  (Sampling) 
#> Chain 4 Iteration: 3390 / 4000 [ 84%]  (Sampling) 
#> Chain 1 Iteration: 3360 / 4000 [ 84%]  (Sampling) 
#> Chain 3 Iteration: 3320 / 4000 [ 83%]  (Sampling) 
#> Chain 2 Iteration: 3550 / 4000 [ 88%]  (Sampling) 
#> Chain 4 Iteration: 3400 / 4000 [ 85%]  (Sampling) 
#> Chain 1 Iteration: 3370 / 4000 [ 84%]  (Sampling) 
#> Chain 3 Iteration: 3330 / 4000 [ 83%]  (Sampling) 
#> Chain 2 Iteration: 3560 / 4000 [ 89%]  (Sampling) 
#> Chain 4 Iteration: 3410 / 4000 [ 85%]  (Sampling) 
#> Chain 1 Iteration: 3380 / 4000 [ 84%]  (Sampling) 
#> Chain 3 Iteration: 3340 / 4000 [ 83%]  (Sampling) 
#> Chain 2 Iteration: 3570 / 4000 [ 89%]  (Sampling) 
#> Chain 4 Iteration: 3420 / 4000 [ 85%]  (Sampling) 
#> Chain 1 Iteration: 3390 / 4000 [ 84%]  (Sampling) 
#> Chain 3 Iteration: 3350 / 4000 [ 83%]  (Sampling) 
#> Chain 2 Iteration: 3580 / 4000 [ 89%]  (Sampling) 
#> Chain 4 Iteration: 3430 / 4000 [ 85%]  (Sampling) 
#> Chain 1 Iteration: 3400 / 4000 [ 85%]  (Sampling) 
#> Chain 3 Iteration: 3360 / 4000 [ 84%]  (Sampling) 
#> Chain 2 Iteration: 3590 / 4000 [ 89%]  (Sampling) 
#> Chain 4 Iteration: 3440 / 4000 [ 86%]  (Sampling) 
#> Chain 1 Iteration: 3410 / 4000 [ 85%]  (Sampling) 
#> Chain 3 Iteration: 3370 / 4000 [ 84%]  (Sampling) 
#> Chain 2 Iteration: 3600 / 4000 [ 90%]  (Sampling) 
#> Chain 4 Iteration: 3450 / 4000 [ 86%]  (Sampling) 
#> Chain 1 Iteration: 3420 / 4000 [ 85%]  (Sampling) 
#> Chain 2 Iteration: 3610 / 4000 [ 90%]  (Sampling) 
#> Chain 3 Iteration: 3380 / 4000 [ 84%]  (Sampling) 
#> Chain 4 Iteration: 3460 / 4000 [ 86%]  (Sampling) 
#> Chain 1 Iteration: 3430 / 4000 [ 85%]  (Sampling) 
#> Chain 2 Iteration: 3620 / 4000 [ 90%]  (Sampling) 
#> Chain 3 Iteration: 3390 / 4000 [ 84%]  (Sampling) 
#> Chain 4 Iteration: 3470 / 4000 [ 86%]  (Sampling) 
#> Chain 1 Iteration: 3440 / 4000 [ 86%]  (Sampling) 
#> Chain 2 Iteration: 3630 / 4000 [ 90%]  (Sampling) 
#> Chain 3 Iteration: 3400 / 4000 [ 85%]  (Sampling) 
#> Chain 4 Iteration: 3480 / 4000 [ 87%]  (Sampling) 
#> Chain 1 Iteration: 3450 / 4000 [ 86%]  (Sampling) 
#> Chain 2 Iteration: 3640 / 4000 [ 91%]  (Sampling) 
#> Chain 3 Iteration: 3410 / 4000 [ 85%]  (Sampling) 
#> Chain 4 Iteration: 3490 / 4000 [ 87%]  (Sampling) 
#> Chain 1 Iteration: 3460 / 4000 [ 86%]  (Sampling) 
#> Chain 2 Iteration: 3650 / 4000 [ 91%]  (Sampling) 
#> Chain 3 Iteration: 3420 / 4000 [ 85%]  (Sampling) 
#> Chain 4 Iteration: 3500 / 4000 [ 87%]  (Sampling) 
#> Chain 1 Iteration: 3470 / 4000 [ 86%]  (Sampling) 
#> Chain 2 Iteration: 3660 / 4000 [ 91%]  (Sampling) 
#> Chain 3 Iteration: 3430 / 4000 [ 85%]  (Sampling) 
#> Chain 4 Iteration: 3510 / 4000 [ 87%]  (Sampling) 
#> Chain 1 Iteration: 3480 / 4000 [ 87%]  (Sampling) 
#> Chain 2 Iteration: 3670 / 4000 [ 91%]  (Sampling) 
#> Chain 3 Iteration: 3440 / 4000 [ 86%]  (Sampling) 
#> Chain 4 Iteration: 3520 / 4000 [ 88%]  (Sampling) 
#> Chain 1 Iteration: 3490 / 4000 [ 87%]  (Sampling) 
#> Chain 2 Iteration: 3680 / 4000 [ 92%]  (Sampling) 
#> Chain 3 Iteration: 3450 / 4000 [ 86%]  (Sampling) 
#> Chain 4 Iteration: 3530 / 4000 [ 88%]  (Sampling) 
#> Chain 1 Iteration: 3500 / 4000 [ 87%]  (Sampling) 
#> Chain 2 Iteration: 3690 / 4000 [ 92%]  (Sampling) 
#> Chain 3 Iteration: 3460 / 4000 [ 86%]  (Sampling) 
#> Chain 4 Iteration: 3540 / 4000 [ 88%]  (Sampling) 
#> Chain 1 Iteration: 3510 / 4000 [ 87%]  (Sampling) 
#> Chain 2 Iteration: 3700 / 4000 [ 92%]  (Sampling) 
#> Chain 3 Iteration: 3470 / 4000 [ 86%]  (Sampling) 
#> Chain 4 Iteration: 3550 / 4000 [ 88%]  (Sampling) 
#> Chain 1 Iteration: 3520 / 4000 [ 88%]  (Sampling) 
#> Chain 2 Iteration: 3710 / 4000 [ 92%]  (Sampling) 
#> Chain 3 Iteration: 3480 / 4000 [ 87%]  (Sampling) 
#> Chain 4 Iteration: 3560 / 4000 [ 89%]  (Sampling) 
#> Chain 1 Iteration: 3530 / 4000 [ 88%]  (Sampling) 
#> Chain 2 Iteration: 3720 / 4000 [ 93%]  (Sampling) 
#> Chain 3 Iteration: 3490 / 4000 [ 87%]  (Sampling) 
#> Chain 4 Iteration: 3570 / 4000 [ 89%]  (Sampling) 
#> Chain 1 Iteration: 3540 / 4000 [ 88%]  (Sampling) 
#> Chain 2 Iteration: 3730 / 4000 [ 93%]  (Sampling) 
#> Chain 3 Iteration: 3500 / 4000 [ 87%]  (Sampling) 
#> Chain 4 Iteration: 3580 / 4000 [ 89%]  (Sampling) 
#> Chain 1 Iteration: 3550 / 4000 [ 88%]  (Sampling) 
#> Chain 2 Iteration: 3740 / 4000 [ 93%]  (Sampling) 
#> Chain 3 Iteration: 3510 / 4000 [ 87%]  (Sampling) 
#> Chain 4 Iteration: 3590 / 4000 [ 89%]  (Sampling) 
#> Chain 1 Iteration: 3560 / 4000 [ 89%]  (Sampling) 
#> Chain 2 Iteration: 3750 / 4000 [ 93%]  (Sampling) 
#> Chain 3 Iteration: 3520 / 4000 [ 88%]  (Sampling) 
#> Chain 4 Iteration: 3600 / 4000 [ 90%]  (Sampling) 
#> Chain 1 Iteration: 3570 / 4000 [ 89%]  (Sampling) 
#> Chain 2 Iteration: 3760 / 4000 [ 94%]  (Sampling) 
#> Chain 3 Iteration: 3530 / 4000 [ 88%]  (Sampling) 
#> Chain 4 Iteration: 3610 / 4000 [ 90%]  (Sampling) 
#> Chain 1 Iteration: 3580 / 4000 [ 89%]  (Sampling) 
#> Chain 2 Iteration: 3770 / 4000 [ 94%]  (Sampling) 
#> Chain 3 Iteration: 3540 / 4000 [ 88%]  (Sampling) 
#> Chain 4 Iteration: 3620 / 4000 [ 90%]  (Sampling) 
#> Chain 1 Iteration: 3590 / 4000 [ 89%]  (Sampling) 
#> Chain 2 Iteration: 3780 / 4000 [ 94%]  (Sampling) 
#> Chain 3 Iteration: 3550 / 4000 [ 88%]  (Sampling) 
#> Chain 4 Iteration: 3630 / 4000 [ 90%]  (Sampling) 
#> Chain 1 Iteration: 3600 / 4000 [ 90%]  (Sampling) 
#> Chain 2 Iteration: 3790 / 4000 [ 94%]  (Sampling) 
#> Chain 3 Iteration: 3560 / 4000 [ 89%]  (Sampling) 
#> Chain 4 Iteration: 3640 / 4000 [ 91%]  (Sampling) 
#> Chain 1 Iteration: 3610 / 4000 [ 90%]  (Sampling) 
#> Chain 2 Iteration: 3800 / 4000 [ 95%]  (Sampling) 
#> Chain 3 Iteration: 3570 / 4000 [ 89%]  (Sampling) 
#> Chain 4 Iteration: 3650 / 4000 [ 91%]  (Sampling) 
#> Chain 1 Iteration: 3620 / 4000 [ 90%]  (Sampling) 
#> Chain 2 Iteration: 3810 / 4000 [ 95%]  (Sampling) 
#> Chain 3 Iteration: 3580 / 4000 [ 89%]  (Sampling) 
#> Chain 4 Iteration: 3660 / 4000 [ 91%]  (Sampling) 
#> Chain 1 Iteration: 3630 / 4000 [ 90%]  (Sampling) 
#> Chain 2 Iteration: 3820 / 4000 [ 95%]  (Sampling) 
#> Chain 3 Iteration: 3590 / 4000 [ 89%]  (Sampling) 
#> Chain 4 Iteration: 3670 / 4000 [ 91%]  (Sampling) 
#> Chain 1 Iteration: 3640 / 4000 [ 91%]  (Sampling) 
#> Chain 2 Iteration: 3830 / 4000 [ 95%]  (Sampling) 
#> Chain 3 Iteration: 3600 / 4000 [ 90%]  (Sampling) 
#> Chain 4 Iteration: 3680 / 4000 [ 92%]  (Sampling) 
#> Chain 1 Iteration: 3650 / 4000 [ 91%]  (Sampling) 
#> Chain 2 Iteration: 3840 / 4000 [ 96%]  (Sampling) 
#> Chain 3 Iteration: 3610 / 4000 [ 90%]  (Sampling) 
#> Chain 4 Iteration: 3690 / 4000 [ 92%]  (Sampling) 
#> Chain 1 Iteration: 3660 / 4000 [ 91%]  (Sampling) 
#> Chain 2 Iteration: 3850 / 4000 [ 96%]  (Sampling) 
#> Chain 3 Iteration: 3620 / 4000 [ 90%]  (Sampling) 
#> Chain 4 Iteration: 3700 / 4000 [ 92%]  (Sampling) 
#> Chain 1 Iteration: 3670 / 4000 [ 91%]  (Sampling) 
#> Chain 2 Iteration: 3860 / 4000 [ 96%]  (Sampling) 
#> Chain 3 Iteration: 3630 / 4000 [ 90%]  (Sampling) 
#> Chain 4 Iteration: 3710 / 4000 [ 92%]  (Sampling) 
#> Chain 1 Iteration: 3680 / 4000 [ 92%]  (Sampling) 
#> Chain 2 Iteration: 3870 / 4000 [ 96%]  (Sampling) 
#> Chain 3 Iteration: 3640 / 4000 [ 91%]  (Sampling) 
#> Chain 4 Iteration: 3720 / 4000 [ 93%]  (Sampling) 
#> Chain 1 Iteration: 3690 / 4000 [ 92%]  (Sampling) 
#> Chain 2 Iteration: 3880 / 4000 [ 97%]  (Sampling) 
#> Chain 3 Iteration: 3650 / 4000 [ 91%]  (Sampling) 
#> Chain 4 Iteration: 3730 / 4000 [ 93%]  (Sampling) 
#> Chain 1 Iteration: 3700 / 4000 [ 92%]  (Sampling) 
#> Chain 2 Iteration: 3890 / 4000 [ 97%]  (Sampling) 
#> Chain 3 Iteration: 3660 / 4000 [ 91%]  (Sampling) 
#> Chain 4 Iteration: 3740 / 4000 [ 93%]  (Sampling) 
#> Chain 1 Iteration: 3710 / 4000 [ 92%]  (Sampling) 
#> Chain 2 Iteration: 3900 / 4000 [ 97%]  (Sampling) 
#> Chain 3 Iteration: 3670 / 4000 [ 91%]  (Sampling) 
#> Chain 4 Iteration: 3750 / 4000 [ 93%]  (Sampling) 
#> Chain 1 Iteration: 3720 / 4000 [ 93%]  (Sampling) 
#> Chain 2 Iteration: 3910 / 4000 [ 97%]  (Sampling) 
#> Chain 3 Iteration: 3680 / 4000 [ 92%]  (Sampling) 
#> Chain 4 Iteration: 3760 / 4000 [ 94%]  (Sampling) 
#> Chain 1 Iteration: 3730 / 4000 [ 93%]  (Sampling) 
#> Chain 2 Iteration: 3920 / 4000 [ 98%]  (Sampling) 
#> Chain 3 Iteration: 3690 / 4000 [ 92%]  (Sampling) 
#> Chain 4 Iteration: 3770 / 4000 [ 94%]  (Sampling) 
#> Chain 1 Iteration: 3740 / 4000 [ 93%]  (Sampling) 
#> Chain 2 Iteration: 3930 / 4000 [ 98%]  (Sampling) 
#> Chain 3 Iteration: 3700 / 4000 [ 92%]  (Sampling) 
#> Chain 4 Iteration: 3780 / 4000 [ 94%]  (Sampling) 
#> Chain 1 Iteration: 3750 / 4000 [ 93%]  (Sampling) 
#> Chain 2 Iteration: 3940 / 4000 [ 98%]  (Sampling) 
#> Chain 3 Iteration: 3710 / 4000 [ 92%]  (Sampling) 
#> Chain 4 Iteration: 3790 / 4000 [ 94%]  (Sampling) 
#> Chain 1 Iteration: 3760 / 4000 [ 94%]  (Sampling) 
#> Chain 2 Iteration: 3950 / 4000 [ 98%]  (Sampling) 
#> Chain 3 Iteration: 3720 / 4000 [ 93%]  (Sampling) 
#> Chain 4 Iteration: 3800 / 4000 [ 95%]  (Sampling) 
#> Chain 1 Iteration: 3770 / 4000 [ 94%]  (Sampling) 
#> Chain 2 Iteration: 3960 / 4000 [ 99%]  (Sampling) 
#> Chain 3 Iteration: 3730 / 4000 [ 93%]  (Sampling) 
#> Chain 4 Iteration: 3810 / 4000 [ 95%]  (Sampling) 
#> Chain 1 Iteration: 3780 / 4000 [ 94%]  (Sampling) 
#> Chain 2 Iteration: 3970 / 4000 [ 99%]  (Sampling) 
#> Chain 3 Iteration: 3740 / 4000 [ 93%]  (Sampling) 
#> Chain 4 Iteration: 3820 / 4000 [ 95%]  (Sampling) 
#> Chain 1 Iteration: 3790 / 4000 [ 94%]  (Sampling) 
#> Chain 2 Iteration: 3980 / 4000 [ 99%]  (Sampling) 
#> Chain 3 Iteration: 3750 / 4000 [ 93%]  (Sampling) 
#> Chain 4 Iteration: 3830 / 4000 [ 95%]  (Sampling) 
#> Chain 1 Iteration: 3800 / 4000 [ 95%]  (Sampling) 
#> Chain 2 Iteration: 3990 / 4000 [ 99%]  (Sampling) 
#> Chain 3 Iteration: 3760 / 4000 [ 94%]  (Sampling) 
#> Chain 4 Iteration: 3840 / 4000 [ 96%]  (Sampling) 
#> Chain 1 Iteration: 3810 / 4000 [ 95%]  (Sampling) 
#> Chain 2 Iteration: 4000 / 4000 [100%]  (Sampling) 
#> Chain 3 Iteration: 3770 / 4000 [ 94%]  (Sampling) 
#> Chain 2 finished in 4792.3 seconds.
#> Chain 4 Iteration: 3850 / 4000 [ 96%]  (Sampling) 
#> Chain 1 Iteration: 3820 / 4000 [ 95%]  (Sampling) 
#> Chain 3 Iteration: 3780 / 4000 [ 94%]  (Sampling) 
#> Chain 4 Iteration: 3860 / 4000 [ 96%]  (Sampling) 
#> Chain 1 Iteration: 3830 / 4000 [ 95%]  (Sampling) 
#> Chain 3 Iteration: 3790 / 4000 [ 94%]  (Sampling) 
#> Chain 4 Iteration: 3870 / 4000 [ 96%]  (Sampling) 
#> Chain 1 Iteration: 3840 / 4000 [ 96%]  (Sampling) 
#> Chain 3 Iteration: 3800 / 4000 [ 95%]  (Sampling) 
#> Chain 4 Iteration: 3880 / 4000 [ 97%]  (Sampling) 
#> Chain 1 Iteration: 3850 / 4000 [ 96%]  (Sampling) 
#> Chain 3 Iteration: 3810 / 4000 [ 95%]  (Sampling) 
#> Chain 4 Iteration: 3890 / 4000 [ 97%]  (Sampling) 
#> Chain 1 Iteration: 3860 / 4000 [ 96%]  (Sampling) 
#> Chain 3 Iteration: 3820 / 4000 [ 95%]  (Sampling) 
#> Chain 4 Iteration: 3900 / 4000 [ 97%]  (Sampling) 
#> Chain 1 Iteration: 3870 / 4000 [ 96%]  (Sampling) 
#> Chain 3 Iteration: 3830 / 4000 [ 95%]  (Sampling) 
#> Chain 4 Iteration: 3910 / 4000 [ 97%]  (Sampling) 
#> Chain 1 Iteration: 3880 / 4000 [ 97%]  (Sampling) 
#> Chain 3 Iteration: 3840 / 4000 [ 96%]  (Sampling) 
#> Chain 4 Iteration: 3920 / 4000 [ 98%]  (Sampling) 
#> Chain 1 Iteration: 3890 / 4000 [ 97%]  (Sampling) 
#> Chain 3 Iteration: 3850 / 4000 [ 96%]  (Sampling) 
#> Chain 4 Iteration: 3930 / 4000 [ 98%]  (Sampling) 
#> Chain 1 Iteration: 3900 / 4000 [ 97%]  (Sampling) 
#> Chain 3 Iteration: 3860 / 4000 [ 96%]  (Sampling) 
#> Chain 4 Iteration: 3940 / 4000 [ 98%]  (Sampling) 
#> Chain 1 Iteration: 3910 / 4000 [ 97%]  (Sampling) 
#> Chain 3 Iteration: 3870 / 4000 [ 96%]  (Sampling) 
#> Chain 4 Iteration: 3950 / 4000 [ 98%]  (Sampling) 
#> Chain 1 Iteration: 3920 / 4000 [ 98%]  (Sampling) 
#> Chain 3 Iteration: 3880 / 4000 [ 97%]  (Sampling) 
#> Chain 4 Iteration: 3960 / 4000 [ 99%]  (Sampling) 
#> Chain 1 Iteration: 3930 / 4000 [ 98%]  (Sampling) 
#> Chain 3 Iteration: 3890 / 4000 [ 97%]  (Sampling) 
#> Chain 4 Iteration: 3970 / 4000 [ 99%]  (Sampling) 
#> Chain 1 Iteration: 3940 / 4000 [ 98%]  (Sampling) 
#> Chain 3 Iteration: 3900 / 4000 [ 97%]  (Sampling) 
#> Chain 4 Iteration: 3980 / 4000 [ 99%]  (Sampling) 
#> Chain 1 Iteration: 3950 / 4000 [ 98%]  (Sampling) 
#> Chain 3 Iteration: 3910 / 4000 [ 97%]  (Sampling) 
#> Chain 4 Iteration: 3990 / 4000 [ 99%]  (Sampling) 
#> Chain 1 Iteration: 3960 / 4000 [ 99%]  (Sampling) 
#> Chain 3 Iteration: 3920 / 4000 [ 98%]  (Sampling) 
#> Chain 4 Iteration: 4000 / 4000 [100%]  (Sampling) 
#> Chain 4 finished in 4912.1 seconds.
#> Chain 1 Iteration: 3970 / 4000 [ 99%]  (Sampling) 
#> Chain 3 Iteration: 3930 / 4000 [ 98%]  (Sampling) 
#> Chain 1 Iteration: 3980 / 4000 [ 99%]  (Sampling) 
#> Chain 3 Iteration: 3940 / 4000 [ 98%]  (Sampling) 
#> Chain 1 Iteration: 3990 / 4000 [ 99%]  (Sampling) 
#> Chain 3 Iteration: 3950 / 4000 [ 98%]  (Sampling) 
#> Chain 1 Iteration: 4000 / 4000 [100%]  (Sampling) 
#> Chain 1 finished in 4929.6 seconds.
#> Chain 3 Iteration: 3960 / 4000 [ 99%]  (Sampling) 
#> Chain 3 Iteration: 3970 / 4000 [ 99%]  (Sampling) 
#> Chain 3 Iteration: 3980 / 4000 [ 99%]  (Sampling) 
#> Chain 3 Iteration: 3990 / 4000 [ 99%]  (Sampling) 
#> Chain 3 Iteration: 4000 / 4000 [100%]  (Sampling) 
#> Chain 3 finished in 4948.7 seconds.
#> 
#> All 4 chains finished successfully.
#> Mean chain execution time: 4895.7 seconds.
#> Total execution time: 4949.1 seconds.
saveRDS(out, "~/simu_rbf.rds")
out <- readRDS("~/simu_rbf.rds")
print(out$converged)
#> [1] TRUE

Inferred dynamics & Rt

source("Figures/figure1.R")
pdf("Figures/figure1.pdf",6,6)
plot(fig1)
dev.off()
#> png 
#>   2
fig1

Pairs Plot

source("Figures/figure2.R")
#> Warning: The `size` argument of `element_line()` is deprecated as of ggplot2 3.4.0.
#> ℹ Please use the `linewidth` argument instead.
#> Warning: `aes_()` was deprecated in ggplot2 3.0.0.
#> ℹ Please use tidy evaluation idioms with `aes()`
pdf("Figures/figure2.pdf",6,6)
plot(fig2)
#> Warning: The dot-dot notation (`..density..`) was deprecated in ggplot2 3.4.0.
#> ℹ Please use `after_stat(density)` instead.
dev.off()
#> png 
#>   2
fig2

Traces & ppcheck

s1 <- plot_traces(out)
pdf("Figures/s1.pdf",10,10)
plot(s1)
dev.off()
#> png 
#>   2
s1

s2a <- plot_ppcheck_At(out,1)
s2b <- plot_ppcheck_At(out,2)
pdf("Figures/s2.pdf",6,6)
plot(s2a)
plot(s2b)
dev.off()
#> png 
#>   2
s2a

s2b

Hyppeparameter pairs

s3 <- plot_hyperpar_pairs(out)
pdf("Figures/s3.pdf",6,6)
plot(s3)
dev.off()
#> png 
#>   2
s3